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The AdS/CFT correspondence has long been used as a tool for understanding 
non-perturbative phenomena in gauge theories because it is an example of a 
'strong-weak' duality: when one side of the duality is weakly coupled, the other 
is strongly coupled and vice-versa. Hence strongly coupled phenomena can be 
studied by looking at the weakly coupled side of the duality. In its original form 
the correspondence proposes a duality between type IIB superstring theory on 
AdS^ x S* 5 and an J\f = 4 supersymmetric Yang-Mills theory in four dimensions. 
In this thesis we investigate proposed duals to QCD itself. Duals to QCD fall into 
two categories: 'top-down' and 'bottom- up'. We take inspiration from both by 
truncating a consistent solution to the type IIB supergravity equations of motion 
(top-down). This model demonstrates dynamical chiral symmetry breaking, has 
a running coupling and contains a holographic description of the vector meson 
sector. By artificially extending the existing U(l) symmetry to SU(2) (bottom- 
up) we then obtain a holographic description of the axial vector sector. We show 
that this model reproduces the masses and decay constants of the lightest mesons 
to the 10% level. By regulating the UV with a sharp cut-off we can reproduce 
the p meson masses to within 2%. Finally we demonstrate that this model can 
be used to reproduce a very good agreement with hadronization data for particle 
production over a range of four orders of magnitude. 
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Chapter 1 
Introduction 



1.1 Introduction Overview 

The advance of physics in the last few centuries has been relentless. One 
of the main themes of this advance has been the desire for unification, the 
ultimate aim being to have one theory that describes all the fundamental 
forces. A major step towards this goal was achieved in the nineteenth cen- 
tury, when separate descriptions of electricity and magnetism were unified 
into a single theory of electromagnetism. It wasn't until the 1930s that 
it was fully appreciated that electromagnetism was an example of a gauge 
theory. 

A gauge theory, in its most general sense, is a model with an invariance 
under a local symmetry of some of the variables in the theory [4,5]. In the 
case of QED, the phases on all the fields can be changed locally, and so long 
as we introduce a gauge field which connects the points of local relabelling, 
the physics remains unchanged. In QED, for example, the quanta of these 
gauge fields are called photons. 
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QED is an example of an Abelian gauge theory: the order in which two 
gauge transformations, GiiGi, are performed is irrelevant. It is possible to 
construct gauge theories that are non-Abelian, meaning that the order in 
which two (or more) gauge transformations are performed does matter. In 
fact it was discovered in the 1970s that this was exactly what was needed 
to describe both the weak force and the strong force, which make up two 
more of the four fundamental forces. The force which we have yet to come 
to is gravity. 

Currently gravity is a bit of an enigma in the particle physics world: 
classically it is well described by the theory of general relativity [6,7], and 
there are no experimental measurements which contradict its predictions. 
However it is hard to test gravity at a single particle level because it is so 
weak in comparison to the other three forces. However almost everyone be- 
lieves that general relativity needs modification; as it stands at the moment 
it cannot be quantized, which would make it unique amongst all the other 
forces. How to quantize gravity has been the bane of quantum physicists 
for generations, and there are currently two prongs of attack: loop quantum 
gravity [8], and string theory [9-11], on which we will spend some further 
time in section fOl However the theme of this thesis is QCD, so it is prudent 
to review some of the features of this non-Abelian gauge theory. 

1.2 QCD 

In the late 1940s and early 1950s, only a few 'elementary particles' were 
known: the proton, neutron, electron, neutrino and photon. Almost the 
entire observable universe consists of just these particles. However, some 
puzzling unstable particles had been seen in cosmic rays. Furthermore, 
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physicists were keen to learn about the nuclear force which was presumed 
to bind protons and neutrons together to form nuclei. This led to the 
construction of larger and larger particle colliders throughout the decades, 
the most recent being the Large Hadron Collider in Switzerland. These 
accelerators have revealed the existence of hundreds of new particles, and 
almost all of them can be classed as hadrons. The discovery of these 
hadrons was not greeted with glee: it seemed incredulous to believe that 
such a large number of particles could all be fundamental. Fortunately, the 
initial skepticism has been validated. A lot of work in the 1950s and 1960s 
has shown that all hadrons (which include the protons and neutrons) are 
composite particles: they themselves are made up of even smaller particles, 
called quarks. It is currently believed that there are six flavours of quark, 
each of which exist in one of three 'colours': red, green or blue. Almost all 
known hadrons can be accounted for by combining a suitable combination of 
quarks in varying levels of excitation. The entire spectrum can be explained 
by allowing gluons (the quanta of the QCD gauge field) to form particles 
too. 

Why quarks are not seen as free entities of their own is a very deep 
question, and is essentially the question of colour confinement. Physicists 
hypothesise that the force between two quarks does not diminish as they 
are separated (contrast this with QED and, say, a positron and electron). 
Therefore it would take an infinite amount of energy to separate two quarks. 
Hence quarks are forever bound into colourless combinations, such as red- 
green-blue or red-antired. These colourless combinations are the hadrons 
we see in particle detectors and cosmic rays. Confinement, although widely 
believed, is yet to be analytically proven. 

When quarks were first discussed in theoretical papers it was not certain 
if they were just convenient mathematical fictions, or whether they were 
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physical entities. An excellent paper [12] predicted what would be seen in 
deep inelastic scattering of electrons and protons if the quarks did physically 
exist. Essentially, the idea was to replicate the idea of Rutherford Scattering 
by firing electrons at the protons and neutrons in nuclei and analyzing 
the resulting scattering pattern. The experiment was performed at the 
Stanford Linear Accelerator Center (SLAC) in 1969. The results matched 
the theoretical predictions, and since then quarks have been embraced as 
physical entities. [13] contains a good review of this experimental evidence. 

1.2.1 The QCD Lagrangian 

The QCD Lagrangian is given by [14] 



C = 4,( ldu Y + gA^t a -m)^--F« v F^ a 
F; v = d,Al-d v Al + gt bc AlAl (1.1) 

Mathematically ^ is a Dirac fermion in the fundamental representation of 
SU(3), A is the gauge field, which is in the adjoint representation of SU(3), 
jabc are structure constants of SU(3), g is a number and m is a mass. 
Physically, ip represents the quarks, with m being their mass. A represents 
the gluons, and g is the QCD coupling constant. Greek letters label space- 
time indices, and Roman letters label gauge group indices. 

The Lagrangian fully describes QCD. It is deceptively simple, and like 
a fractal that on close inspection is ever more complex, so too is QCD. For 
example, it is not immediately clear from f 1 1.1 1) how to calculate hadron 
masses, decay constants or scattering cross-sections. In what follows we 
choose to mention those aspects which will be of relevance in this thesis. 
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1.2.2 Renormalization 

The tricky subject of renormalization is dealt with in many textbooks [14,15] 
and is the result of efforts made in the 1940s when it was realised that simple 
calculations, even in QED, resulted in divergent integrals. It is now accepted 
that these divergences, if finite in numbecl, can be reabsorbed into the pa- 
rameters of the Lagrangian. Hence quantities such as the quark mass, m, 
the QCD coupling constant, g, and the field strengths, A^, as they appear 
in the Lagrangian in equation (11. ip are infinite. These quantities are called 
the bare values, and when we speak of the quark mass, QCD coupling con- 
stant etc. we are actually referring to the redefined renormalized quantities. 



To properly construct a systematic method of renormalization, we have 
to define a renormalization scale which is the energy at which we define 
our renormalized quantities (m, g, A^). We then see what happens to these 
quantities as we vary the energy scale at which we are working. It turns 
out that the quantities vary as we change the energy scale. What must be 
remembered at all times is that the physics is independent of our choice 
of renormalization scale: experimentally measurable quantities are com- 
pletely blind to which renormalization scheme the theorist has chosen. The 
implication of this statement is explored in the following section, via the 
renormalization group equation. 

1.2.3 The Renormalization Group Equation 

The first step towards calculating hadron masses, decay constants and cross- 
sections from (11.11) is the calculation of correlation functions. In fact this is 

*A theory is, by definition, rcnormalizable, if there are only a finite number of diver- 
gences. If there are an infinite number, the theory is said to be non-renormalizable. For 
details on how we know whether a theory is renormalizable or not, see [14] 
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the most difficult step. 

Let us consider the bare connected n-point correlation function, in 
4> 4 theory! It is defined as 



r ( g\x 1 ,x 2 ,...,x n ) = (0\T(j)o(x 1 )(f)o(x2)...<j)o(x n )\0) 



;i.2) 



where <po{x) refers to the bare fields, before renormalization. They are 



A/2 



(x). The 



related to the renormalized fields by a rescaling: <po{x) - 
bare quantity, is a function of the bare coupling constant, go, and some 

(n) 

cutoff, A. The renormalized quantity, T R is a function of the renormalized 
coupling constant g and the renormalization scale M. 



T^(k l ;g,M) = Z { ; /2) T%\k l ;g ,A) 



(1.3) 



The bare theory is independent of the renormalization scale, M, so we can 
write 

( M ^)(^" n/2)r(n)(fc ^' M) ) = ° (L4) 
It is easy to show that this equation can be rewritten as 



> , 9 _ d n 



r^(k f ,g,M) = 



;i.5) 



with the definitions 




1 



M 



M 



dg 
dM 
d\nZ 4 

dM 



(1.6) 
(1.7) 



theory is one of the simplest examples of an interacting quantum field theory. 
Its Lagrangian is C — ^(d^cj)) 2 — ^m 2 (j) 2 — jifi 4 - More details can be found in [14]. 
However, its exact nature is unimportant here. We use it because its renormalization 
group equation is particularly simple. 



7 



A little thought will also convince us that (3 and 7 can only be functions 
of g: by dimensional analysis they cannot depend upon M, and are clearly 
independent of either n or {xi}. Using the method of characteristics, we 
can solve (11.51) : 

T^\k i; g,M) = T^\k i; g(p);M)ex V ^ J%(g(p))^ (1.8) 

with the characteristic equation p^jj^- = f3(g(p)), and the initial condition 
g(p = 1) = 9- 

Equivalent equations to (11.61) and (11.71) can be derived for all quantum 
field theories, and the definitions of (3 and 7 are identical to those definitions 
in (USD and (Ol . 

Next we turn to the significance of (3 and 7. 

1.2.4 The beta function 

We defined the beta function, /3(g) in (11.61) . and from its definition we can 
see that the beta function tells us how the renormalised coupling constant 
changes as we vary the renormalisation scale. The strength of the coupling 
constant at any particular energy is an important quantity: it determines 
when perturbation theory is valid. Perturbation theory is the most suc- 
cessful tool we have for solving quantum field theories: it was used very 
successfully in QED, so it would be excellent if we could use the same tech- 
nique for QCD. 

At low energies QCD is not weakly coupled: the coupling constant is 
large and perturbation theory is not applicable. But what about at high 
energies - perhaps we could use perturbation theory then? 

To see what happens, let us assume that at high energies the theory 
is weakly coupled, and check that this doesn't lead us to a reductio ad 
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absurdum argument. 

If g is small, we can approximate the beta function at these energies by 
a Taylor series expansion: 



(3(g) =A + Bg + Cf + Df... (1.9) 

where g is the deviation of g from its value at p = M. We define the value 
of g at p = M as g R : g R = g(p = M). 

Any perturbative deviation from g R will involve some sort of interaction. 
With each interaction, at least one power of the coupling constant will enter. 
Hence we can set A = 0. 

B,C, D . . . can be calculated quantitatively^; the method is detailed 
in [14]. Here we are only concerned with the qualitative effects of the the 
leading coefficient of (11.91) . The first non-zero contribution to (II. 9p is canon- 
ically called f3\. There are three possibilities: 

1. f3t > 

2. f3t = 

3. (3i < 

For (3i(g) > 0, the coupling constant becomes small in the infrared region 
of the theory. This is the case with QED, and allows us to successfully use 
perturbation theory at every-day energies. However, in the high energy, 
ultraviolet regime, the coupling constant gets larger and larger, eventually 
reaching the point where perturbation theory is no longer valid. Fortunately, 
in QED, the energy level where this occurs is much much higher than any 
reachable today. Furthermore, new physics is expected to appear a long 
time before we reach such energies. 



tin QCD with 6 flavours of quark, B = C = and D = -jf^ [14] 



For (3\{g) = 0, the coupling constant doesn't run, and so the bare cou- 
pling constant is equal to the renormalised coupling constant for all energies. 
This is the case for many supersymmetric theories. It is also the case with 
scale invariant field theories and conformal field theories: if the 
coupling constant doesn't run, and there is no other source of scale in the 
theory (e.g. a mass scale could be introduced by having massive quarks), 
then the theory is defined to be scale invariant. A conformal field, however, 
is invariant to metric rescalings (ds 2 — > X 2 ds 2 for arbitrary, real A). All con- 
formal field theories are scale invariant (the metric can be used to measure 
distance, in a conformal field theory we can always rescale the metric, so 
it is impossible to unambiguously define a distance), however not all scale 
invariant theories are conformal. Known examples are rare though [16-19]. 

For Pi(g) < 0, the coupling constant gets smaller at large energies, and 
conversely gets large at small energies. This is the case for QCD (see next 
section). So although we cannot use perturbation theory at low energies for 
QCD, we can use it for high energies; typically those found at large particle 
accelerators in Switzerland, Germany and California. 

The beta function for QCD 

In 1973, the one-loop contribution to (3(g) for QCD was calculated to be 

with ri/ being the number of quark flavours, and N c deriving from the 
SU(N C ) gauge group. There are widely believed to be six flavours of quark, 
so /3i(g s ) < 0, meaning that at high energies the quarks interact weakly. 
They are said to be asymptotically free. 
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Figure 1.1: Plot of the solution to the first loop beta function for QCD 
We can go further, and solve f ll. 10D : 

9s( P ) 2 = (1.H) 

whose solution is plotted in figure 11.11 Note that p enters the solution, 
for it is the only dimensionful parameter available to make the log term 
dimensionless. So in fact we can swap the M derivative with a p derivative 
in the definition of the (3 function in fll. 6ft : 



! U 77T7 771 77' — ' ^ = TTi L12 

oM oiogM aiogp 



This reparameterization will come in useful in the next section. 

Figure 11.11 brings to our attention two important features of QCD. 
Firstly, g gets smaller and smaller as p — > oo. So QCD is indeed asymptot- 
ically free, and we can use perturbation theory at high energies. 

Secondly, figure 11.11 demonstrates the anomalous breaking of scale 
invariance. If we were to look at the massless QCD Lagrangian, ie. (11.11) 
with m set to zero, we would see that there is no dimensionful parameter. 
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We say that, classically, massless QCD is scale invariant. However, this 
would be a puzzle: phenomenologically there is clearly an inherent mass 
scale in QCD: the hadronic masses are not continuous, so any one could be 
used to define a scale. On studying the (3 function we would be reassured. 
Figure [T7TI clearly introduces a scale into the theory, commonly called Aqcd- 
It is normally defined to be the energy at which the coupling constant 
becomes infinite. In figure 11.11 this is approximately 246 MeV (although a 
more careful analysis, including higher order corrections to /3(g) sets Aqcd 
at ~ 215 MeV). For all energies where g s > 1, perturbation theory is no 
longer valid. Quantum effects have broken the classical symmetry: we say 
that the scale invariance has been broken anomalously. 

1.2.5 The 't Hooft Expansion 

Equation (11.101) provides an ideal opportunity to talk about the 't Hooft ex- 
pansion, a particularly clever way of performing perturbation theory about 
the strong scale Aqcd- We have just shown that at and around Aqcd> the 
coupling constant is greater than one, which means a perturbative expan- 
sion such as in (11.91) would be invalid. However in [20,21] it was pointed 
out that there is a second dimensionless parameter in QCD: that of 'iV c ' in 
the gauge group SU(N C ). 't Hooft showed that gauge theories simplify at 
large N and that they have a perturbative expansion in terms of 1/N. We 
summarise his analysis, and that of [22], in this section. 

We stated in (ll.lOj) that the /3-function for a pure SU (N c ) gauge theory 
is given by 

Pigs) = ~y n °Tw + ° {9l) (L13) 

so the leading terms remain the same if we let N — > oo, so long as we keep 
A = g 2 s N fixed (one can show the higher order terms also stay the same in 
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this limit). This limit is known as the 't Hooft limit. 

The 't Hooft expansion was originally formulated in terms of a U(N C ) 
gauge theory, with all matter in the adjoint representation (although this 
can be generalised to include fundamental matter). Let us assume, that 
as in QCD, all three-point terms are proportional to g s and all four-point 
terms are proportional to g 2 , so the Lagrangian takes the schematic form 

£ ~ Tr(d$A) + g.j jk Tr($& j Q k ) + g 2 a d ikl Tr(<i> i $ j $ k $ l ) (1.14) 

for some constants c tJ and d l ^ kl . We can rescale the fields by $j = g s $i, to 
reach 
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Tr{d$id$i) + c ijk Tr(^^ k ) + d^Tr^^i) (1.15) 



where \/g\ M = N/X. Now we want to know what happens in the limit 
N — > oo, A = fixed. Although 1/g 2 tends to infinity, this is countered by 
the fact that the number of components in the fields also goes to infinity. 

A U(N C ) theory with matter in the adjoint representation can be repre- 
sented as a direct product of a fundamental and an ant i- fundamental field 



$*, as in figure fL2l where we have drawn two contributions to the vacuum 
amplitude. 

What is the power of N and A associated with each diagram? From 
(11.15!) we can see each vertex will have a coefficient proportional to 
N/X, and each propagator will have a coefficient proportional to X/N. In 
addition, each closed loop will a factor of N c since we have to sum over all 
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Figure 1.2: Two examples of contributions to the vacuum amplitude in dou- 
ble line notation. The diagram on the left is planar and hence contributes 
N 2 . The diagram on the right is of genus one and hence contributes N° 



indices in the loop. Hence we have 



N 

Vertex : (1.16) 
A 

Propagator : — (1.17) 
Loop : iV c (1.18) 



So a diagram with V vertices, E propagators and F loops will include a 
factor 

N N-E + F X E-V = N X X E-V (Llg) 

where x is the Euler characteristic of the diagram. It is a topological 
invariant, depending only on the genus (number of holes), g, of the surface: 

X = 2 - 2g (1.20) 

Therefore the perturbative expansion of any diagram in field theory can be 
written as 

oo oo oo 

£a^E c ^ = £^ 2 - 2 ^w (i.2i) 

9=0 i=0 g=0 

where f g is some polynomial in A. In the large iV limit, we see that any 
computation will be dominated by those diagrams which are topologically 
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equivalent to a sphere or plane (ie. no holes, genus = 0). 

If we identify 1/N as some generic coupling constant, g p , we can see 
that (11.211) looks like a perturbative expansion about g p . The resemblance 
between (11.2ip and perturbative string theory is one of the strongest moti- 
vations for believing that string theories and field theories are related [22], 
in particular by a 1/N expansion. 

In particular, the AdS/CFT correspondence, which is the inspiration for 
this thesis, and which we turn our attention to in section [2j is formulated 
at large N. It is hoped that in the future 1/N corrections will make the 
correspondence even more applicable to £77(3) quantum chromodynamics. 

1.2.6 Anomalous Dimensions 

In section 11.2.41 we looked at (3(g) when g was small. Now we turn our 
attention to g > 1, the strongly coupled regime. We can no longer calculate 
(5(g) explicitly, but we can consider the qualitative possibilities. 

Of most interest is when (3(g) is either positive or negative in the weakly 
coupled regime, but higher order corrections mean that (3(g) has non-trivial 
zeros: the two possibilities are shown in figure 11.31 




Figure 1.3: Possible forms of the (3 function with non-trivial zeros: (a) 
ultraviolet-stable fixed point; (b) infrared-stable fixed point 

A beta function of the form 11.3( a) is weakly coupled at low energies. 
At higher energies the coupling constant grows, but only upto g*. Once it 
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reaches the coupling constant stays constant. This is a potential, but 
unproven resolution of the Landau pokt! in QED: the Standard Model has 
such a pole at 10 34 MeV [23]. However this analysis assumes a monotonic 
beta function. If the beta function of QED has a non-trivial zero, then the 
Landau pole would not exist. 

Explicitly, near the fixed point of figure fLBT a). the (3 function can be 
approximated by 

P~-B(g-g*) (1.22) 



or 



which has the solution 



g~-B(g-g.) (1.23) 



d\ogp 



m = 9* + c(jY (1.24) 



This has important implications for the exact solution of the renormalization 
group equation (11.81) . For sufficiently large p, the integral in the exponential 
of fll.Sp will be dominated by those values of p where g(p) is close to g*. 
Then 



G(p) « G (^)exp [-(log ^2(1- 7 G?*))] (1.25) 
/ 1 \ 1 ~ 7 ( 9 *) 

- C {-2) (!-26) 



So the two-point function returns to a simple scaling law, as might be ex- 
pected in a naive dimensional analysis argument. But there is an important 
difference: we would expect the power law to be p~ 2 , but instead it is 



§A Landau pole means that the coupling constant becomes infinite at finite energy. 
Technically QCD has a Landau pole at ~ Aqcd, but the phrase is usually used to refer 
to non-asymptotically free theories only 
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p-2(i-7(s*)) _ The complex interactions of the quantum field theory have af- 
fected the law of rescaling. Finally, since the fixed point is reached at large 
p, it is known as an ultraviolet fixed point. 

An analysis of graph 11.3( b) follows similar lines. Near the fixed point 
the P function can be approximated by 



(3 ~ +B{g - g.) 



1.27) 



which has the solution 



g(p) = g* + c ( 



p_ 

M 



B 



:i.28i 



This means that for sufficiently small p, the integral in the exponential of 
(11.81) will be dominated by those values of p where g{p) is close to g*, giving 



-(bg ^2(l-7(a,)) 



C 



P~ 



1-7(9.) 



;i.29) 
;i.3o) 



Because this occurs at small p, this fixed point is called an infrared fixed 
point. 

Once again, the complex interactions of the field theory have affected 
the law of rescaling. For this reason, the 7(g) function is commonly known 
as the anomalous dimension, even if there is no fixed point in the theory. 



1.2.7 Mass gap 

Another feature of QCD, and yet to be proved rigorously by anybody, is 
that of a mass gap. A quantum field theory is said to have a mass gap 
if the energy spectrum has a positive greatest lower bound, but does not 
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include zero. This is thought to be very closely linked to the property of 
confinement, which is simply a statement that gluons cannot exist on their 
own, but only in colourless bound states. Lattice gauge theories have shown 
to the satisfaction of most that quarkless QCD exhibits such a phenomena, 
but solving this problem with due mathematical rigour still remains one 
of the Millennium Prize Problems. Experimentally it is certainly true that 
there is a mass gap in QCD. 

1.2.8 Chiral Symmetry Breaking 



We have already seen in section 11.2.41 that massless QCD anomalously 
breaks its classical scale invariance. There is a second classical symme- 
try to the massless QCD Lagrangian: that of chiral symmetry, meaning 
that the left and right handed quarks transform independently. 

The Lagrangian of massless QCD is given by taking equation ( 11.11) and 
simply setting m = 0: 



C = 4, {id vl v + gA a ^t a ) $ ~ \F a , u F^ a 



;i.3i) 



fll.3ip possesses a chiral symmetry which can be seen if we write ip as 
and choose the chiral representation of the gamma matrices: 



7 



v 1 V 



r 



\ 



-a 1 



/ 



then (suppressing the gauge field which doesn't affect the analysis), equation 
fll.3ip becomes 



£ = (d - a. y) ipL + (d + a. y) ^ 



R 



;i.32) 



18 



In this form it is clear that we can perform separate global flavour transfor- 
mations on ipi and iprt.- Since the quarks are in an SU(2) isospin multiplet, 
we write these transformations as SU{2)l <g> U(1)l and SU(2)r <g> U{1)r 
respectively. The Lagrangian is unchanged by these transformations: it is 
chirally invariant. (The analysis still holds if we turn on the gauge field A^.) 

However, this is a symmetry of the Lagrangian which is broken sponta- 
neously by the choice of vacuum. The most obvious manifestation of this 
is the absence of a parity-doubled spectrum. If the chiral symmetry was 
respected, there would be a positive parity hadron for every negative parity 
hadron. This is not what is seen in nature. For example, the proton has 
isospin ~, spin | and a positive parity. Its mass is 938 MeV. Its parity part- 
ner, catchingly called the N(1535) Sn, has a mass of 1535 MeV [24]. We 
conclude that chiral symmetry is broken at a quantum level, and the easi- 
est way to account for this is by assuming that QCD spontaneously forms 
a quark condensate. In other words, the QCD vacuum gives a non-zero 
vacuum expectation value to the scalar operator 



What is the physical interpretation of this VEV? The up and down quarks 
are very light, and so it costs little energy to create a quark-antiquark pair. 
However, the binding energy released by a bound quark-antiquark pair is 
large. All (11.331) is saying is that the process of forming a quark-antiquark 
pair and then binding them together is exothermic: we get out more energy 
then we put in. 

The vacuum expectation value fll .33j) means that we can no longer 
perform independent gauge transformations on the left and right handed 



(0\QlQr + QrQl\0) ^0 




19 



spinors: the only flavour transformation we still can make is one where we 
simultaneously make the same transformation on both and We say 
that the SU(2) l x SU {2)r flavour symmetry has been spontaneously broken 
to SU(2) V . 

Goldstone's theorem [25] states that whenever a continuous symmetry 
of a quantum field theory is spontaneously broken, a massless particle will 
appear. In breaking SU (2) L x SU (2) R to SU (2)y, we started with six group 
generators, and ended up with three: according to Goldstone's theorem we 
would expect three massless particles to appear in the spectrum of QCD. 

However, the Lagrangian of (ll.3ip is not exactly that of real QCD. The 
proper Lagrangian of QCD has massive quarks and is given in ( II. ip . Hence 
QCD does not have an exact SU{2) L x SU{2) R symmetry QCD. But be- 
cause the up and down quarks are almost massless, there is an approximate 
symmetry, so a perturbative approach about m q pa may still be useful. 

Looking at the hadron spectrum there are three suspiciously light hadrons, 
namely the pions, whose masses are about one fifth of the next lightest 
hadron. Furthermore they have the correct parity to be created by the 
axial isospin current, j^ 5a : 

The coupling between the pions and the vector axial current is defined as 

(Ap)\Ql,l 5 T a Q(0)\0) = -tf nP ,6 ab (1.35) 

where T a is a generator of the broken symmetry group and f n is a number 
with dimensions of mass. It is called the pion decay constant. For an SU(2) 
isospin symmetry, meaning two flavours of quark, we expect three types of 
pion. 
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Figure 1.4: The vector current two-point function, IIy(— q 2 ). The grey blob 
represents all possible one particle irreducible diagrams that can be inserted. 



f n in equation (11.351) is a measure of the extent of the chiral symmetry 
breaking. If there is no chiral symmetry breaking its value is zero. In QCD 
it can be determined from measuring the decay rate of 7r° particles [24]. It 
is found to be 92 MeV. 

Meson decay constants 

It is not just the pion that has a decay constant associated with it. All 
pseudoscalar mesons, P, have a decay constant Fp, defined by 

(0\A li (0)\P{q))=iF P q ll (1.36) 

And the vector meson decay constantly are defined as: 

(0|Q(0) 7 M Q(0)|V>; A)) = iFymyel (1.37) 

where p and A are the momentum and polarization state of the vector meson 
V(p; A). e M is the corresponding polarization vector. 

Experimentally, Fp can be measured from leptonic decays of the ap- 
propriate meson [24]. Fy can be measured from the decays of tau lep- 

^vector mesons have two decay constants: the transverse polarization decay constant, 
Fy, and the longitudinal polarization decay constant Fy. Only the longitudinal decay 
constant can be determined experimentally. 



21 




1 



111' 



2 




V 



Figure 1.5: One contribution to the vector current two-point function 
tons [26-28]. 

If we represent the vector current two-point function as figure 11.41 then 
one contribution to the vector two-point current will be that of figure 11.51 

It can be shown [29, 30] that at large N, the mesons form an infinite, 
stable spectrum and that the vector two-point current is given by the sum 
over all meson and glueball resonances (with the correct quantum numbers), 
such as that in figure 11.51 



The axial vector current two-point function at large N, Tlp(— q 2 ), is almost 
identical: 



Axial Current Anomaly 

fll.3ip also possesses a U(1)l <8> U {1)r symmetry which we have yet to men- 
tion. Whereas the 577(2) symmetries mixed the up and down quarks, the 
£7(1) symmetries act to add a phase to the entire spinor. Classically, we have 




(1.38) 




(1.39) 
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Figure 1.6: One of the diagrams that leads to the axial vector anomaly in 
QCD 

independent U(l) symmetries for both the left and right handed symmetries. 
These phases always cancel out since the fermions are always bilinear. The 
presence of the non-zero VEV (11.331) also breaks this classical U(1)l<8>U(1)r 
symmetry to leave a single U(l)y symmetry. And once again, Goldstone's 
theorem applies: a spontaneously broken symmetry implies the presence of 
a massless boson. Here the broken symmetry has one only one generator, so 
we'd expect one massless boson. A look at the known QCD mass spectrum 
shows that the only suitable candidate, in terms of quantum numbers, is the 
if (958) meson. However it has a disappointingly high mass. Even allowing 
for the small masses of the quarks, we'd expect our boson to have a mass 
similar to that of the pions. 

The explanation of this dichotomy puzzled physicists for a long while, 
and wasn't explained until 1986 [31]. It is quite easy to show [14,32] that 
due to the necessity of regularization and renormalization (in particular 
triangle diagrams such as figure ll.6p the classical axial symmetry of QCD 
(<9/xJ m5 = 0) is broken by quantum effects to 

W = -4r^ VF ^ ( L4 °) 

1D7T 

The obvious question to ask is what is the source of such an anomaly, and 
the answer lies in the complicated field of instantons [31-34]. There are 
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an infinite number of distinct QCD vacua, whom differ from one another 
by their topological equivalence class: 7T 3 (S 3 ) = 7L. Furthermore, it can be 
shown that whenever a correlation function tunnels in and back out of a 
vacuum in a different topological class, it picks up a contribution to the 
axial symmetry: 



where Q is defined as the topological charge, and can be shown to be equal 
to the number of left handed fermions minus the number of right handed 
fermions: Q = Ul — Ur. An instanton is defined as those solutions for which 



Probably the simplest way to think of instantons is as isolated sources of 
axial symmetry breaking, and they are responsible for making the if (958) 
meson heavy. More detailed analysis can be found in the sources quoted 
earlier. 



Now that we've investigated some of the complex features of QCD, a pru- 
dent question would be to ask whether we have the theoretical tools required 
to calculate experimental quantities such as hadron masses and decay con- 
stants. These are very difficult quantities to determine. Current approaches 
can be divided into four classifications. 

• Perturbative QCD 

Feynman diagrams are used highly successfully in QED to calculate 
cross-sections. The correlation function is expanded in terms of a 
small parameter, a ~ 1/137, and then each contribution is calculated 
order-by-order. Beyond one loop the calculations can get very in- 




(1.41) 



Q= 1. 



1.3 Solving QCD 
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tractable, and it requires a great deal of effort and care. Fortunately, 
a is very small, and so higher order contributions are also small. The 
equivalent expansion in QCD is a s ~ 0.118 at 91.2 GeV [35], which 
means that the expansion does not converge as rapidly. In addition, 
as we discussed earlier, we can only use perturbative QCD at high 
energies, so for many interactions, it is not valid. 

• Lattice QCD [36, 37] 

The best established approach to non-perturbative QCD is the use 
of huge supercomputers. Continuous spacetime is represented by a 
discrete four dimensional lattice, with lattice spacings of a. The QCD 
Lagrangian (II. ip is then reformulated for a discrete spacetime, and 
the calculations are made for as small a lattice spacing as possible. 
For the final answer, the limit a — > is taken. Lattice QCD takes a 
lot of human time and effort, and it is not always clear what happens 
when the limit a — > is taken. 

• Effective Theories 

Mirroring the techniques of perturbation theory, theories are writ- 
ten down which mimic QCD for some limited parameter space. And 
within this parameter space, there will be a parameter which can be 
used as a perturbation parameter. Examples are chiral perturbation 
theory [38-40], which uses the light quark masses (at low momentum) 
as perturbation parameters, and heavy quark effective theory [39,41], 
which uses the inverse of a heavy quark mass as a perturbation pa- 
rameter. These approaches can be very useful, but are limited by the 
extent of parameter space to which their model is valid. 
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• 1 /N expansion 

Once again, mirroring the techniques of perturbation theory, experi- 
mental quantities are calculated as a perturbative expansion in terms 
of the N in the SU(N) of the gauge group of QCD, in line with the 
argument presented in section 11.2.51 In theory this is applicable to all 
of QCD, but the perturbation parameter, (1/3) 2 , is not very small, 
so the series does not converge as quickly as one might like. More 
importantly, it is not clear how to go about calculating the 1/N cor- 
rections. Hence there is little knowledge of the error in a tree-level 
calculation. The AdS-CFT correspondence falls under this classifica- 
tion. Before we investigate it we need to go over some basic aspects 
of string theory. 

1.4 String Theory 

The AdS / CFT correspondence is discussed in detail in chapter [21 It de- 
scribes a duality between a gauge theory and a string theory, so we take 
this opportunity to review some of those aspects of string theory relevant 
to this thesis. 

1.4.1 Premise 

String theory starts by assuming that fundamental particles are not point- 
like, but have a finite length. In other words, they are little strings. 

We reach the Euler-Lagrange equations of motion for point-like particles 
by minimizing their world-lines. By analogy, we assume that the equations 
of motion for strings will be found by minimizing their world areas. The 
original formula for the area of the sheet, known as the Nambu-Goto action 
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[9,10], is 



S = T I dadr^X 2 X' 2 - {X.X') 2 



(1.42) 



where 



X^ 



dX^(a,r) 



do 



(1.43) 



The X^ are the bosonic fields in the spacetime. (11.421) does not contain 
fermionic fields; this is addressed later. 

The square root in (11.421) makes the mathematics of optimization awk- 
ward. Instead we can introduce a new variable, h a p which will be the metric 
on the string world sheet. This leads to the conventional 



Since the derivatives of h do not appear in (jl.44p . it is not a real field, and 
so can be eliminated, leading back to ( 11.421) if required. 



(I1.44p is not satisfactory to describe the real world. Firstly, we need 
fermions in the theory, and secondly it turns out that (ll.44p contains a 
tachyon. However it is possible to supersymmetrize (jl.44j) [9]: 



S[X, 4] = -— d 2 o{rf h gilv d a X»d h X») + a^^vf^^f) (1.45) 



And with a lot of work [9] it is possible to show that so long as X M and ^ 
exist in ten dimensions, (11.451) contains no tachyons or anomalies. 

1.4.2 World-sheet Boundary Conditions 

There are several types of string: their properties depend on the nature of 
the boundary conditions of the fields on the world-sheet. For the moment we 




(1.44) 
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Type 


Details 


I 


All strings are charged under an SO (32) gauge 

cjA/m m pt rv 

o y iiiinv v l y 


IIA 


Chiral, meaning that parity is a good symmetry. 

l\Fn crancrp wmmptn/ pifnpr f npnrAy odIa/ pnntairm 

1 i \J &jCX Oj Xlllllv^ Ll V v3 1 uXlv^X . OU (. 1±V7 V v 1_ V Will y IjUllt dilllo 

gravity 


IIB 


Non-chiral, both in the fermionic sector and the 

era n (tp ^pptor 


HO 


Hptprotic mpanin£ r that thp left-movprs arp 
bosonic, and the right-movers are fermionic. Has 
an SO (32) gauge symmetry 


HE 


Heterotic, with an E 8 x E 8 gauge symmetry 



Table 1.1: The five consistent types of superstring theory 

only consider closed strings (open strings enter when we consider D-branes 
in section 11.4.31) . We can choose to make our strings orientable or non- 
orientable. If a string is orientable, we can tell which way we are travelling 
along a string. If it is non-orientable, we can't. There are further choices 
we can make with respect to the fermionic boundary conditions, and in fact 
going through every possibility took researchers a lot of time. In the end, 
there are five consistent superstring theories, summarised in table ll.lt 

What is the spectrum of closed string theory? This is a detailed and 
complex subject that is visited by many sources (see [9,11,42-45] for exam- 
ple) so here we just briefly summarise the results. 

The spectrum contains the metric, G^, the scalar dilaton $, and a two 
index antisymmetric tensor B^ u : 

• The metric needs little explanation, and plays the usual role as the 
device for measuring distance and angles. 

• The dilaton is a very important field in string theory: it generically 
appears in the action as e*, and measures the Euler characteristic 
(the number of holes and handles) of the string worldsheet. Hence 
the quantity e* plays the role of the theory's coupling. When open 
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and closed string sectors are combined the Yang-Mills coupling from 
the open string sector has g\ M = e*. 

• B, known as the Kalb-Ramond or NS-NS B- field, is the generaliza- 
tion of the one-form potential of electrodynamics. For point-like 
particles moving in an electromagnetic potential the action is given 
by J A^dx^. For closed strings moving in a Kalb-Ramond potential, 
the action is given by J B^dx^dx 1 ' . 

All of the fields listed above are NS-NS fields (Neveu-Schwarz-Neveu- 
Schwarz). R-R (Ramond-Ramond), NS-R and R-NS fields also exist. The 
NS-R and R-NS fields are not a concern to us - they are simply the fermionic 
partners of the bosonic sector. However the R-R fields are sourced by D- 
branes [9] and so are also present in the spectrum. They enter as field 
strengths of antisymmetric fields as shown in (11. 46ft and (11.471) . 

For the AdS/CFT correspondence applied to 3 + 1-dimensional field 
theories, ten-dimensional type IIB string theory is of central importance, 
and in particular its low-energy limit where strings become point-like and 
string theory becomes supergravity. There exists no completely satisfactory 
action for the type IIB supergravity, since it involves an antisymmetric field 
C4 with self-dual field strength F 5 . However, it is possible to write an 
action involving both dualities of C4, and then impose the self-duality as 
a supplementary field equation. In this way one obtains (see for example 



[42-44]) 




1 



(1.46) 
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where the field strengths are defined by 



F x = dC, H 3 = dB, F 3 = dC 2 , F 5 = dC 4 , 
F 3 = F 3 - CH 3 , F 5 = F 5 - X -A 2 A H 3 + l -B A F 3 , 



(1.47) 



and we have the additional self-duality condition *F 5 = F 5 . 

However, for the AdS/CFT correspondence, there are only three fields 
switched on: the graviton, the dilaton, and the five- form field strength F 5 , 
so we can disregard most of the fields in (11.461) and instead use: 



action from the string frame to the Einstein frame, giving an action in the 
usual Einstein-Hilbert form [46]. It is in this frame that the stress-energy 
tensor has its usual meaning. 

There is more to the string theory story than that indicated by table fTTTl 
As well as the 1+1 dimensional strings, there are higher dimensionful objects 
too, called D-branes. Their importance was first recognised in the mid- 
1990s, and they form an important part of the AdS-CFT correspondence. 
We look at them next. 

1.4.3 D-branes 

Upto this point we have been conducting a purely mathematical exercise. 
It was hoped that string theory would provide a unique theory of the four 
fundamental forces of nature. However, we have already stated that there 
are five consistent versions of superstring theory. This is disillusioning, and 
was the point at which the string community had reached in the early 1990s. 
Fortunately, it was shown that all five theories were linked by a series of 




(1.48) 



where we have redefined the metric (y/Ge 



yG) to transform the 
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simple dualities, meaning that there was only one theory after all, called 
M-theory. One of those dualities is T-duality. 

T-duality is a symmetry between small and large distances. The 'T' 
stands for topology, since T-duality transformations work on spaces in which 
at least one dimension has the topology of a circle. 

Let us consider a string that is wrapped around the circle, ie. it is closed. 
It's energy spectrum is given by 

9 n ui K ,„ 
m> ~ ^ + ^ d-49) 

The first term is the Kaluza-Klein energy, and the second term is the wind- 
ing number. (11.491) is invariant under the following transformation: 

R->a'/R, n-^uj (1.50) 

Therefore a theory with only closed modes will be invariant under (ll.SOp . 
But what if a theory contained open strings too? These aren't invariant 
under f 1 1 . 5 j) . The open strings will have a Kaluza-Klein contribution to 
their energy spectrum, but no winding term. Hence T-duality will change 
the spectrum of the open strings. 

Under T-duality, the boundary conditions of the open string change: 

d a2 X»(a u (0, J)) = <— X»(<t u (0, 0) = C (1.51) 

We go from Neumann boundary conditions (LHS) to Dirichlet boundary 
conditions (RHS). For a long time this was result was ignored by the string 
community because the Dirichlet boundary condition is not Lorentz invari- 
ant - a special surface in the 10D spacetime has been picked out. It was 
Polchinski who first embraced this result, and extolled the modern explana- 
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tion of this effect. The surface c M is indeed special; it represents the surface 
of another object, hitherto undiscovered in string theory. These objects 
are called D-branes, the 'D' standing for 'Dirichlet', and the 'brane' as in 
membrane, or surface. By the T-duality, all open strings must start and 
finish on a D-brane. 

Just as closed strings can be considered fluctuations of the background 
geometry, open strings can be considered as fluctuations of the D-branes. 
The D-branes are solitonic objects, and thus exist in the low energy limit 
of string theory - supergravity. 

1.4.4 D-brane action 

It was shown [47] in 1989 that the action 



for some arbitrary D-brane gives the equations of motion for all the back- 
ground fields associated with any D-brane. This form of an action was 
first written down by Born and Infeld when considering a particular non- 
linear generalization of electromagnetism [48], and was expanded upon by 
Dirac [49]. All this occurred a long time before the advent of string theory 
and D-branes, but for historical reasons (11.521) is still called the Dirac-Born- 
Infeld (DBI) action. Let us briefly explain each field in turn: 

• $ is the usual dilaton. 

• G' is the induced metric on the brane's surface. The precise definition 
of this concept is given in section 13.2.11 but it is the equivalent of 
the yj— detg term in the standard Einstein-Hilbert action for general 




(1.52) 



relativity. 
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• B' is the induced Kalb-Ramond field on the brane's surface. The 
precise co-ordinate definition can be found in [9], but all it represents 
is the two-form field which has been induced on the surface of the 
D-brane as a result of the Kalb-Ramond field in the bulk. 

• F is the usual electromagnetic field potential for a single D-brane, and 
is attributable to the open string sector. In fact the B' fields and F 
fields interact, and it can be shown [50] that only the combination 
F + B' is gauge invariant. 

In many publications, it has been conventional to re-define F as 2na'F, 
and for clarity we follow this convention in this thesis. Furthermore, the 
Kalb-Ramond field will be switched off for all the examples in this thesis. 
Hence the DBI action for the rest of this thesis is taken to mean 



There will be cases though when we turn off the electromagnetic field tensor 
F in (TL53D too. 

For completeness, we add that it is possible to add a Chern-Simons term 
to flL52l [51,52]: 



However this is a topological term that is only relevant if we are considering 
vacuums with non-trivial topologies. In all that follows, we implicitly as- 
sume that the vacuum is topologically trivial, and hence we will not include 
a Chern-Simons term in any subsequent action. 




(1.53) 




J M 



/2 v / det(G + B> + F) + i / Tr exp(27ra'F 2 + B 2 ) A ^ C q 



(1.54) 
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Chapter 2 
AdS-CFT 



The discovery of the AdS/CFT correspondence was built upon a whole 
raft of previous discoveries and conjectures. Here we ignore the historical 
steps that led to its discovery, and simply introduce it, as if by magic, as it 
stands today. On one side we have a type IIB string theory in the geometry 
induced by an infinite stack of D3 branes. The other side is a conformal 
gauge theory with an infinite number of colours. First we address the D3- 
brane construction. 

2.1 String description 

Our starting point is to put N D3-branes in the centre of a ten dimensional 
space, which is otherwise empty. There will be three types of interaction: 
those amongst the open strings on the branes; those amongst the closed 
strings in the bulk; and those between the open and closed strings. We can 
write the effective action as 

S = •S'bulk + 'S'branc + ^int (2-1) 
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Now we consider the limit in which we send the string length to zero; l s —>■ 
(a' — > 0). All other dimensionless parameters (the string coupling constant, 
g s , and N) we keep fixed. Doing this means the coupling of the strings to 
each other (~ g s ot' 2 ) goes to zero, and hence the strings do not interact: 
Si n t is turned off, and we are left with two decoupled systems: classical ten 
dimensional gravity in the bulk, Sbuik, an d a four dimensional gauge theory 
on the surface of the branes, SWane- 

What is the exact nature of the four dimensional gauge theory? The ends 
of the open strings are labelled by the branes they are attached to. The 
strings are orientated, so there are two types of strings stretching between 
any two branes: one whose left end is on the first brane and whose right end 
is on the second brane, and vice versa. These labels are called Chan-Paton 
factors, and in the case of a coincidental stack of branes, as we have here, 
there is a clear re-labelling symmetry of the Chan-Paton factors. 

There are N 2 possible string relabellings, which fills the adjoint repre- 
sentation of the U(N) Lie group. The strings are free to move around on 
the surfaces of the branes, so the relabelling can be performed at any point: 
hence the U(N) symmetry is a local one. In other words it is a gauge theory, 
and the Chan-Paton labels can be interpreted as labelling the gauge charge 
of each end of the string. 

In the large N limit, the propagators of adjoint fields in a U(N) and 
SU(N) gauge theory differ by a vanishing term, proportional to iV -1 [22], 
so when N — > oo, we can treat U(N) and SU(N) gauge theories as identical. 
The spectrum of strings is given by the dimensional reduction of an N = 1 
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;auge multiplet in ten dimensions to four dimensions: 




(2.2) 



D=4,Af=4 



where the left hand side of (12.21) is defined in ten dimensions with one 
supersymmetry generator (J\f = 1), and the right hand side is defined in 
four dimensions, with four supersymmetry generators (A/" = 4). The are 
gauge fields, \1/ are gauginos, and A represent complex scalars. Six of the 
ten-dimensional gauge field degrees of freedom become the three complex 
scalars in the four- dimensional reduction. 
The one-loop (3 function for this theory is 
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-N 



-N-N 







(2.3) 



with the first term coming from the gauge fields, the second from the gaugi- 
nos, and the third from the scalar fields. Supersymmetric non-renormalization 
theorems [53] can be used to show that there are no higher order contri- 
butions to the j3 function, and thus it is zero for the full quantum theory. 
A theory which has no massive fields, and whose coupling constant doesn't 
undergo dimensional transmutation, is a conformal field theory (CFT). This 
means the theory is invariant to rescalings of the metric. 

To conclude: on the string side, we have two decoupled systems: ten- 
dimensional supergravity; and an SU(N) conformal gauge theory. 
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2.2 Supergravity description 

Let us consider the system from a different point of view. D-branes are 
massive objects, and hence warp the spacetime that they are put in. ft is 
possible to show that there is a D3 brane solution to supergravity [54] which 
is 

6 

ds 2 = Z- 1 'S ilv dx li dx v + Z 1 ' 2 dxad%a 

o=l 

Z = 1 + ^ R 4 = ^g s l 4 s N [ *F 5 = N (2.4) 
r Js 5 

R characterises the typical curvature scale of the gravitational solution, 

and N is the number of D3-branes in the spacetime. 

Now let's take the same low energy limit as we did before (l s — > 0). From 
the point of view of an observer at infinity, there will be two types of low 
energy excitations. The first type would be any massless excitation in the 
bulk with a suitably large wavelength. The second is a bit more subtle, and 
is a result of g tt not being constant. All observations made at infinity will be 
redshifted, by a factor Z^ 1 ^, when compared to the observations made at 
some fixed position r. So we can have very highly excited states in this low 
energy system, just so long as they are close enough to r = 0, because by the 
time they have reached some observer at infinity the redshift factor will be 
so large that all excitations will appear low energy to that observer. And as 
we tend closer and closer to the limit l s — > 0, the two systems will decouple: 
the excitations near r = because the gravitational well is insurmountable; 
and the supergravity excitations because the wavelengths will be very large 
compared to the gravitational extent of the branes. 

In conclusion, we have supergravity in the bulk, and in the near horizon 
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we have 




where we have been able to go from (12.41) to (12.51) because in the limit we are 
considering, r -C R. Equation (I2.5p is simply the geometry of AdS 5 x S 5 , 
with the radius of the 5-sphere given by R. 

We see that from both the string description, and the supergravity de- 
scription, the system decouples into supergravity in flat space, and some- 
thing else. This suggests the conjecture that the 'something else' are al- 
ternative descriptions of exactly the same thing: this is the essence of the 
AdS/CFT correspondence: 

jV = 4U(N) SYM in 3+1 dimensions 
is dual to type IIB superstring theory on AdS$ x S" 5 

It is sometimes easy to get lost in the mathematical details of what we 
are doing, but it is worth pausing to consider exactly what we are claiming. 
The SYM side of the duality is defined in 3+1 dimensions: the supergravity 
side exists in 9+1 dimensions! In some sense, six of the ten dimensions 
on the supergravity side are superfluous! This type of duality is classed as 
a holographic duality, taking its name from holograms which use a two di- 
mensional object to depict a three dimensional image. The first holographic 
dual was proposed for black-holes, with the postulate that the entropy of a 
black-hole is given by its area. This meant that all the information contained 
within a black hole is encoded in a lower dimensional surface. Similarly we 
are claiming here that all the information in the ten dimensional supergrav- 
ity can be equally well explained by a four dimensional field theory. 
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2.3 Matching the Symmetries 

What symmetries do each side of the description have? On the AdS$ x 5 5 
side we clearly have two symmetry groups: the 50(6) of the 5-sphere, and 
the symmetry of the AdS§ spacetime is SO (4,2) [22]. What does this relate 
to on the super Yang-Mills side? Since it is M = 4 there are four gauginos, 
hence there is an SU(4:)r = 50(6) symmetry. It is also conformal in all 
four dimensions: the conformal group in four dimensions has a symmetry 
group SO (4,2). So we can identify corresponding symmetries on either side 
of the correspondence: 





M = 4 SYM 


AdS 5 x 5 5 


SO(4,2) 


conformal group 


spacetime symmetry 


SU(4)=SO(6) 


R- symmetry 


5 5 isometry 



2.4 Regime of usefulness 

In what regime are our approximations valid? On the supergravity side, the 
two regimes of interest decoupled when the wavelength of the supergravity 
excitations were much larger than the string length: 

1 < 4 ~ 9y M N ~ 9 S N (2.6) 

L s 

But on the field theory side we need the 't Hooft coupling to be much less 
than one if the theory is to be in the perturbative regime: 

1 » 9ymN -9sN^^- (2.7) 

Clearly the two inequalities are incompatible, which leads us to the conclu- 
sion that when the gravity dual is weakly coupled, the field theory must be 
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strongly coupled. We can be bolder still, and make the conjecture that when 
the gravity side is strongly coupled, the field side will be weakly coupled. 
We propose that the duality is a strong-weak duality. 

This is why this duality is so exciting to those studying strongly cou- 
pled gauge theories. As explained earlier, solutions to the non-perturbative 
regime of such theories have evaded capture for decades. This duality gives 
us hope that instead of having to try and solve the nightmare of a strongly 
coupled theory, we can instead re-parameterise the theory in terms of a 
weakly coupled supergravity theory, which we can solve using perturbation 
theory! 

2.5 Energy-radius Duality 

Ignoring for the moment the SU(4)ji = SO(6) part of the duality, we have 
a 4+1 spacetime dual to a 3+1 field theory. Conventionally the 'extra' 
fifth dimension of the spacetime is denoted with r, and here we consider its 
physical interpretation. 

The AdS side has the metric 

ds 2 = r~ 2 dr 2 + r 2 ri llv dx^dx v (2.8) 

But the CFT is invariant under x M — > e a x fl . If the CFT side is invariant 
under this transformation, then so too must be the AdS side. Hence r must 
scale as r — > re~ a . It transforms as an energy. 

This can be understood from the point of view of an observer at infinity 
(the boundary of AdS$ x S 5 ). To her, a photon emitted from the centre of 
the space is red-shifted. If the observer moves in from infinity to some finite 
r, the photon is less red-shifted. So the energy-radius duality indicates that 
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instead of just considering the duality as being constructed on the boundary 
of the AdS$ x S* 5 space, we are free to look at the duality at any value of 
r: each different choice of r simply represents the same field theory at a 
different energy scale. 

From here onwards, the r coordinate is referred to as the energy scale. 
The small r limit is called the infrared (IR), and the high r limit is called 
the UV (ultraviolet). 

2.6 Explicit Correspondence 

The correspondence will only be useful if there is a prescribed way of cal- 
culating physical quantities on each side. The way to do this was first 
suggested in [55], and the postulate is that the generating functional for 
correlation functions with some source 4>o(x) is equal to the supergravity 
partition function whose fields <f>{x) tend to (j)o(x) at the boundary of the 
AdS 5 x S 5 space: 



So, in simple terms, the boundary values for supergravity fields correspond 
to sources for the field theory operators. 

The LHS side of (12. 9ft is independent of the global symmetries of the 
conformal field theory. So, the entity J d A x (p (x)O(x) must be a singlet 
under all the symmetries of the dual field theory. This is how we go about 
identifying which fields in the supergravity match with which operators in 
the CFT, which we address next. 




string [00 ( x )\ 



(2.9) 
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2.7 Field-operator Matching 

We showed in the previous section that an operator O on the conformal 
side corresponds to putting in a field <fi on the gravity side. 

The example canonically presented is to consider the action of a scalar 
field, 0, of mass m on the AdS$ x S 5 side: 

S = J d 4 x dr dtt 5 {rf v d^d v (\) - m 2 2 ) (2.10) 

The simplest solution is obtained by assuming that is a function of r, the 
radial direction of the AdS$ x S* 5 space only. Then the solution is given by 



Ar -(4-A) + Br -A ( 2n > 



with m 2 = A(A - 4j 

In fact in general, a bulk scalar with mass m 2 corresponds to a scalar 

Solving 



it, 



operator on the conformal side with dimension A = — 2± v4 + 
□0 = gives two independent solutions, r A_4 and r _A : 0(r) ~ Ar A ~ 4 + 
Br~ A . We can then identify [55-58] S as the source of the operator O 
and A as the VEV for the same operator. As discussed in the previous 
section, the most natural way [55] to combine these quantities to form a 
dimensionless quantity is J d A x (A X B) = J d 4 x 4>qO. 

We also stated in the last section that J d 4 x(po(x)0(x) must be a singlet 
under all symmetries of the conformal field theory. Furthermore, it must 
have mass dimension since it is part of an exponential, so 4>o(x)0(x) must 
have mass dimension 4. Clearly such a combination of A and B satisfies 
these requirements. 



*This formula can be generalized to m 2 = (A — p)(A + p — 4) in the more general 
case of the action of a p-form [22] 

refers to the full dimension of the operator (i.e. classical + anomalous) 
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Let us consolidate these ideas by looking at the example of the gaugino 
bilinear AA. It is a Lorentz scalar of dimension 3: therefore we need a su- 
pergravity scalar field with mass dimension 1. It is also a ten dimensional 
representation of the SU(4) R-symmetry. We find such a field in the reduc- 
tion of type IIB strings on AdS$ [55] if we set m 2 = — 3(3 in (12.101) . Then 
the solution becomes 

0(r) = Ar^ + Br" 3 (2.12) 

m c , . 

= 7 + ^3 ( 2 - 13 ) 

where we have just redefined the constants A and B. Supergravity fields 
do not scale under four- dimensional conformal transformations, so m must 
have dimension 1 and c must have dimension 3, which corresponds with m 
being a source (i.e. a mass) for A A and c being its VEV (or condensate). 

The careful reader may ask at this point about existence and uniqueness 
when it comes to operator-field matching. Fortunately that question has 
been answered in [59]: for every operator there exists a corresponding field 
which is unique. However, as is the way with existence-uniqueness theories, 
there is no easy way to match operators with fields: each one has to be done 
by hand. And in any gauge theory there are an infinite number of operators. 
When we come to consider duals to QCD, a lot of the work is in deciding 
which operators of QCD, from an infinite choice, are to be included in the 
holographic field theory. 



"'"AdS has negative curvature, and so we can have m 2 < without faster-than-light 
travel 
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2.8 Deformations and Adding Flavour 

So far we have explored the AdS/CFT duality as it was first presented, 
basically as a mathematical construct. However the promise of describing 
a strongly coupled gauge theory in terms of a weakly coupled one was very 
exciting to students of QCD. Perhaps, it was thought, this would be a useful 
way to make calculations in the strongly coupled regime of QCD. However, 
there were some immediate obstacles to overcome: the CFT field theory 
present in the duality had many characteristics which makes it distinctly 
different to QCD. It is conformal, supersymmetric, strongly coupled in the 
UV, and has an infinite number of colours. In addition there aren't any 
quarks. The addition of quarks is addressed in the next chapter. 
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Chapter 3 

Flavouring the AdS$ x S 



All the matter discussed in chapter [5] was in the adjoint representation of the 
SU(N) colour group. This is because both ends of the open strings start and 
end on the same type of brane - a D3-brane - and so are indistinguishable. 
All the Chan-Paton factors are in effect colour indices. Quarks have both 
colour and flavour indices, and are in the fundamental representation of 
the gauge group. So we need strings in our geometry which can take on 
these properties. A lot of work, by a lot of contributers, showed that this is 
achieved by introducing D7 branes [45,60-76]. In this section we summarise 
the salient points of their work. 

3.1 The D7-Brane Probe 

A brane is massive, and will distort the spacetime in which it is placed. At 
the moment this is an unwanted complication for us, so we shall use the 
probe limit, which means that we ignore the gravitational effects of the new 
brane. This approximation, known as the quenched approximation, holds if 
Nf <C N c , where Nf is the number of probe branes, and iV c is the number of 
D3 branes. This approximation is the same as introducing a test charge in 



45 





Xq 


Xi 


%2 
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D3 
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D7 


X 


X 


X 


X 


X 


X 


X 


X 







Table 3.1: Choice of D3 and D7 embedding. Filled directions are marked 
with a cross. Unfilled directions by a dot. 

electrodynamics. Work has been done on unquenched spacetimes [70,77-80], 
but this is not too much of a concern for us here. 
The AdSs x S 5 metric can be written as 

ds 2 = ^ v dx»dx v + ^Y1 dr i ( 3 - X ) 

i=i 

A D7-brane probe can be chosen to fill the x directions and four of the r 
directions, as illustrated in table 13.11 

Introducing the D7 branes has allowed two new types of string to exist 
in this system: 

• Strings stretched between the D3 and D7 branes - these have one 
colour index and one flavour index, and so can be identified as quarks 

• Strings with both ends on a D7 brane - these have two flavour indices, 
and so are identified as mesons 

This is, of course, in addition to the strings with both ends on a D3 
brane ('gluons') and closed strings ('gravitons') which were present before 
the D7 branes, and are still there. 

If the D3 and D7 stack are coincidental in the r$ — plane, there is a 
conformal symmetry in the system. If they are not coincidental, then we 
have introduced a scale into the system, and the conformal symmetry is 
broken; this is shown in figure 13.11 In breaking the conformal symmetry, all 
the fermions and scalars of the M = 2 multiplet acquire a mass. 
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/•D7 




D3 



Figure 3.1: If we place the D7 brane(s) at any non-zero distance from the 
D3 branes, we break the conformal symmetry 

The exact embedding of Nf D7-branes has been chosen so that the 
directions of and are orthogonal to the D7-branes' world-volume (table 
13.11) . There is now an M = 2 chiral supermultiplet connecting the D3 and 
D7 brane stack, and this interacts with the M = 4 supermultiplet on the 
surface of the D3 branes. How to know which symmetries are preserved 
and which are broken in some setup of Dp and Dp' branes is not an obvious 
one, and involves some detailed calculations [81]. In this particular case, the 
M = 4 supermultiplet is broken completely by the presence of the D7 brane, 
leaving just the M = 2 supermultiplet. When the D7 brane is separated 
from the D3 stack, the scalar associated with the r§ and r 6 directions gains 
a VEV and gives a mass to the quarks: 

Q^5gQ — ► rnQQ as A 56 acquires a VEV (3.2) 

However when when we have separated the D3-D7 system, there is still 
a U(l) symmetry because of our freedom to re-parameterise the r$ — r$ 
coordinates. This U(l) symmetry corresponds to the axial symmetry of a 
single-flavour QCD model. As explained in section H. 2. 81 in QCD this sym- 
metry is spontaneously broken by the chiral condensate (qq). We don't see 
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a Goldstone boson though, because of the presence of instantons. However, 
in the 't Hooft limit (TV — ► oo with gyyi-^ fixed), the instanton effect is 
negated, and so the Goldstone boson becomes light once again. 

3.2 Particle spectrum in the probe limit 

The excitations of the D7-brane in the directions perpendicular to its world 
volume correspond to scalar fields, which are singlets under the colour group 
and are in the adjoint of the flavour group: these states include mesons. 

Since the theory is supersymmetric, there will also be other states, in- 
cluding fermions, living on the brane, but these are not explored in this 
thesis. 

The presence of supersymmetry also means that a fermion condensate 
cannot be expected to form (we check this explicitly in the next section). 
Hence chiral symmetry will not be spontaneously broken, as it is in QCD. 
Non-super symmetric deformations of the correspondence can overcome this 
drawback, and in chapter H] we address a particularly important example of 
such a deformation. 

3.2.1 Quark mass and condensate 

In this section we analyse how the probe D7 brane lies when introduced into 
the D3-brane setup. We will show that its asymptotic solution is holographic 
to the quark mass and condensate. 

The full action for a probe D7-brane in the D3-brane setup is 

S D7 = -T 7 J d 8 ^-det (P[G] + 2ira'F) + T 7 J P[C (4) ] A F A F 

(3.3) 

where G is the ten dimensional metric, and F is the field strength of the 
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gauge field U(l) living on the D7 brane. P[G] is the pullback of the ten 
dimensional metric onto the surface of the D7 brane. It is defined as 



P[G] - %%G ab (3.4) 



with a, (3 running over the coordinates on the D7 brane, and a, b running 
over the whole spacetime's coordinates. 

For now, we neglect the U(l) field and the Wess-Zumino term, leaving 
us with 

S D7 = -T 7 J d 8 ^-det (P[G]) (3.5) 

We can neglect the Wess-Zumino term because it is only relevant for gauge 
fields with a vector index on the 5-sphere. In other words they carry a 
supersymmetry charge. Such states are important in the full dual, but for 
this thesis they are unimportant - we are trying to provide holographic 
descriptions of QCD, which does not (to the best of our knowledge) contain 
any supersymmetric states. We will re-introduce the U(l) field in section 
13.2.21 It will turn out to be holographic to the vector meson sector. 
The DBI action for (13.51) is 



' ^P 3 L.f du 5 2 , du 6 2 \ R 2 ( du 5 2 du 6 2 \ 



where the ten-dimensional metric has been written as 



u 2 R 2 
ds 2 = —rj^dx^dx" + — (dp 2 + p 2 dn 2 3 + du\ + du 2 6 ) (3.7) 



with u 2 = p 2 + u\ + tig. The x direction is perpendicular to the D7's world 
volume, so fluctuations in x correspond to the meson spectrum. We come to 
this next, but first we address fluctuations in the p direction - this governs 
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how the probe brane lies within the spacetime. Using the U(l) freedom 
we have to re-parameterise 7/5 and uq, let us set u$ = 0, and look at the 
equation of motion for 7/5. Expanding the square root to second order, the 
Euler-Lagrange equation of motion is 



whose solution tends to m + cp~ 2 as p —> oo. 

From (I3.6p we can see that 7/5 has dimensions of mass. Hence the pa- 
rameter m corresponds to a mass for the qq operator, and the parameter c, 
with mass dimension three, corresponds to a VEV for the qq bilinear. 

Now we can investigate whether this system does or does not allow a 
quark condensate. For the field theory to be consistent, we require two 
things of the D7 flow: 

• The field theory must have a unique description at any given energy: 
hence the sum u\ + p 2 must be monotonic as p varies from to oo 

• The metric (13. 7ft is symmetric under p — > —p, which implies that well- 
behaved solutions must obey t/' 5 (0) = and 7/5(0) 7^ ±00. Were it 
otherwise, the brane would have a kink, which would give an infinite 
contribution to the action 

Figure 13.21 shows three flows for m = 1 with c = —1,0 and 1. As 
expected for a supersymmetric theory, a condensate is not allowed to form. 

If the reader is not satisfied with this numerical explanation, we can 
provide an analytical proof. (13. Sp can be solved analytically, for all p. Its 




(3.8) 
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Figure 3.2: D7-brane flows for m = 1 with several different condensate 
values in AdS^ x S 5 background. Only the c = solution allows a consistent 
theory 



solution is 



i 2/3 B 1 
uJp) = A + - F(6, -(2 + y/3)) x 

£2 + i2 /3 r 2) (£4 _ z 2/ 3jB 2 r 2 + ^4/3 r 4) (3 9) 



with (f) = arccos( ^ 2 "^^ / 3~ 1 + ^^^ ). F(</>, m) is the incomplete elliptic integral 
of the first kind, using the convention defined in Mathematica version 5.2. 

The function is complex unless B = 0, meaning that the only well- 
behaved solution is when w 5 is equal to a constant. This is the solution 
with a zero valued condensate. By allowing m^Owe have clearly broken 
the U(l) chiral symmetry that existed before. We showed in section 11.2.81 
that QCD breaks chiral symmetry both by the spontaneous formation of a 
chiral condensate, and to a lesser extent, by the mass of the quarks. In this 
pure AdS system, the only way we can break chiral symmetry is by making 
the quarks massive. 
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Meson mass spectrum 

How to calculate the meson spectrum in this background has been detailed 
in [62]. We summarise the results here. 

The meson spectrum corresponds to x-dependent fluctuations of the 
fields u 5 and u 6 on the D7-brane world volume. To find the spectrum, 
the mesons are treated as small fluctuations on the brane surface, meaning 
that the DBI action can be expanded to second order, with no interaction 
between the meson fields. Hence we can treat the fields as freely moving 
plane waves, which will allow us to explicitly calculate the spectrum. 

There are two directions in which the fields and Uq can fluctuate: 
perpendicular to the p-dependent flow of the brane, and parallel to the 
p-dependent flow. We use the U(l) symmetry to write this as 

U5 = m + r) u 6 = 0(x,r) (3.10) 

The induced metric on the D7 brane can be written as 



ds 2 = ^^n^dxV + -^jjdp 2 + f^^l (3-11) 



with p 2 = u 2 — L 2 . L is the distance between the centres of the D7 brane 
and the D3 brane stack. It is therefore proportional to the quark mass: the 
constant of proportionality is unimportant in this analysis. It just means 
that we will predict the masses of the mesons relative to one another, rather 
than relative to some absolute scale. 
The DBI action is then given by 

Sdbi = ~T 7 J d 8 £ V^fetG (l + 2(7ra'R) 2 ^j^(d cX d d x + d c <j>d d <f>fj 

(3.12) 
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The equations of motion for <fi and x are identical, so we write the generic 
field as $, and its equation of motion is 



(p 2 + L 2 ) 2 " p 

We assume that the mesons will be non-interacting and plane-wave, so we 
make the ansatz $ = f{p)e lk ' x , with k 2 = —M 2 (M being the mass of the 
meson). We then look for normalisable solutions to ( 13. 13f) : these correspond 
to consistent field theory solutions. 
(13. 131) can be solved exactly by: 

f(p) = (p 2 + L 2 y a 2 F 1 (-a, -a + 1; 2; -p 2 /L 2 ) (3.14) 

2-Fi(a, b; c; z) is the regular hypergeometric function, defined as 

2 F x {a^-c-z)^^f^ Z ^ (3.15) 

where (a) n is the rising factorial, a is as yet undetermined. We can de- 
termine it by considering our requirements on f(p): the solution must be 
normalizable and real. To stop (p 2 + L 2 )~ a 2-^1 (— «, — ol + 1; 2; — p 2 /L 2 ) run- 
ning off to ±00 as p — > 00 we need to terminate the series at order p~ 2a . A 
little maths shows that this can be done by setting 

n = a — 1 (3.16) 

in which case the hypergeometric function terminates at order (p 2 /L 2 ) n . 
Hence ~ p~ 2 as p — > 00. So the final solution to (13.131) is 



f(p) = {p 2 + ^pi " 1. ""J ^ V/^ 2 ) (3.17) 
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The condition ( 13. 161) means that the mass spectrum is discrete: 

2L 

M 2 n = —y/(n + l)(n + 2) (3.18) 

Equation (13.181) shows that in pure AdS^ x S 5 , the meson mass spectrum 
is proportional to the quark mass. This is at odds to what we'd expect in 
QCD: even if the quarks were massless we'd still expect the mesons of QCD 
to be massive. This is another sign that this geometry does not support 
a quark condensate. A meson in QCD gets the majority of its mass from 
interactions with the non-zero quark condensate: when this condensate is 
zero, its mass becomes proportional to the quark mass. 

Numerical Solutions 

We have been lucky in being able to solve (13.1 3D exactly: in the more 
complicated geometries to come, we will not be able to, and so it will be 
helpful to go over the alternative numerical technique, first detailed in [45], 
using an example where we can check that the answers are what we'd expect. 

The task is to find those values of M n in (13.131) for which f n (p) is a 
normalizable, real and smooth solution. 

1. Define the numerical solution f n (p) as a function of M n . To define the 
numerical solution requires us to set boundary conditions, so we make 
the appropriate choice in the UV: for the meson spectrum we require 
f n {oo) = p~ 2 . So now f n is a function of M n and p: f n = f n (p, M) 

2. We now want those solutions which are finite in the IR and smooth 
across the origin: otherwise the solution will have a kink across p = 
0. We could do this by altering M slightly each time, and then re- 
plotting the solution. A better method is to define a new function 
g(M) = sgn [d p f(0, M)\ and plot g(M) over some suitable range. 
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Figure 3.3: Numerical calculation of the meson spectrum in AdS^ x S* 5 . The 
dots are the analytical values. 



3. Where the value of g(M) changes sign indicates a stable flow: an 
example plot for (13.131) with L — 1, M — \[2 is shown in figure 13.31 



Figure 13.31 shows that the numerics agree with the analytical solution. 
This will be important - in future we will not be able to find an analytical 
solution, so we have to be confident that the numerical solutions can be 
trusted. To give an even clearer idea of what is happening around one of 
the normalizable solutions to (13.131) we plot /(r) for three different values of 
M around -v/2 in figure 13.41 It can be seen that as we pass through M = V2 
the solution flips from — oo to +oo, so the only normalizable solution in this 
neighbourhood is at M = v2. 



3.2.2 Vector Meson Spectrum 



We analysed the scalar meson sector in section 13.2.11 Now we re-introduce 
the U(l) gauge field that we neglected in equation (13. 5ft : 



S D7 = -T 7 J d^V-det (P[g] + 2tt«'F) 



(3.19) 
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Figure 3.4: Plot of three solutions around the M = \f2 solution. 

The process is identical to that in section 13.2.11 once we have made the 
following ansatz for the holographic gauge field: 

A» = g(p)e ik - x e» (3.20) 

e M is the polarization vector, defined only in the Minkowski directions. g(p) 
is determined from the equations of motion, just as f(p) was for the scalars. 
In this geometry, the mass spectrum of the vector mesons is exactly the 
same as that for the scalars: 

2L 

M 2 n = —y/(n + l)(n + 2) (3.21) 

This is to be expected since they both form part of an M = 2 hypermultiplet. 
In the more complicated geometries that we will study later the M = 2 
hypermultiplet will be broken, and then the masses will differ. 
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3.3 Breaking the Supersymmetry 

As it stands at the moment, supersymmetry is still exact, and hence a 
fermion bilinear condensate cannot be formed. In other words, chiral sym- 
metry is not broken, which is a necessity for any holographic description 
of QCD. To induce a breaking of the supersymmetry, we must deform the 
gravity side in a non-supersymmetric and non-conformal way. There are 
many examples of such deformations [82-108]. The one of most relevance 
to us is that considered in [109] and we address it next. 
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Chapter 4 



The Constable-Myers 



Geometry 



We showed in chapter [3] that D7 probe branes can be used to include quarks 
in holographic duals. However, the resulting theory still had a massless 
quark condensate. Here we take a step back, and before introducing quarks, 
we consider a gravity dual where supersymmetry has already been broken. 

4.1 The Background 

The Constable- Myers geometry, first studied in [109], is a consistent solution 
to the type IIB supergravity equations of motion. We state the solution 
first, and then describe its properties. In the Einstein frame, the geometry 
is given by: 



w 



w 4 



,4 6 



ds 2 = H~ 1/2 f^dxj + R 2 H X I 2 f( 2 ~ 5 ^ 



(4.1) 



where H = f - 1, / 



w 4 + fr 4 
w 4 -b 4 



(4.2) 



58 



The dilaton and four-form are 




(4.3) 



There are formally two free parameters, R and b, since 



5 = R 4 /2b\ A 2 = 10 - 5 2 



(4.4) 



The metric (14. lj) has a singularity at = 6 4 ; we loosely expect this to 
correspond to the presence of the central stack of D3 branes [69]. 

The dilaton and five-form field strength are a function of r, and in the 
UV limit (r — > oo) the geometry returns to AdSs x S 5 . However, due to the 
dependence on r, the interior of the field theory is strongly deformed. 

The set of solutions are described by two dimensionful parameters, b 
and R. However, by looking at (14.11) we can see that if b — 0, then the 
solution returns to AdSs x S 5 , with radius R. Hence b determines the 
scale of the conformal symmetry breaking. We can identify 6's holographic 
partner explicitly by looking at the form of the dilaton in the UV limit 
(r — »■ oo): 



A has scaling dimension 4 and is a singlet under the 5-sphere isometry. 
Therefore it must correspond to an operator of dimension 4 which is an 
R-singlet. There is only one candidate (as pointed out in section 12.71 we 
are guaranteed to have such a unique candidate). It must be (TrF 2 ). This 
deformation completely breaks the supersymmetry of the gauge theory (ie. 
from M = 4 to M = 0). In addition, a complex VEV does not make 



e 




(4.5) 
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much sense, so b must be greater than or equal to 40 _1 ^ 8 i?. Below this, the 
geometry is not well-defined. 

Meanwhile R determines gy M iV in the field theory, as is usual in the 
correspondence. So R specifies if the field theory is in the strongly coupled 
or weakly coupled regime. 

4.2 Confinement and Glueballs in the Constable- 
Myers Geometry 

Before we introduce quarks into this geometry, it is worth pausing to con- 
sider some of the properties already present in this theory, namely confine- 
ment and the glueball spectrum. 

4.2.1 Confinement 

Confinement means that quarks are consigned to exist in colour-neutral 
combinations. A single quark cannot exist on its own. One way of showing 
this is to calculate the energy of a quark-antiquark as a function of their 
separation. In QED, the energy between a positron-electron pair would 
monotonically decrease to zero as we separated the two. In QCD, we'd 
expect the energy to increase as we tried to separate a quark-antiquark 
pair, showing that it would take an infinite amount of energy to produce a 
solitary quark. 

In [110] it is argued that this effect can be studied holographically by 
looking at the Nambu-Goto action for a fundamental string in the geometry 
dual to the field theory. Specifically, the expectation value of a Wilson loop 
in the field theory is dual to the fundamental string in the holographic 
background: 



60 



L 



L 

Figure 4.1: Wilson loop used to calculate the quark-antiquark force. The 
vertical direction represents Euclidean time. The horizontal direction indi- 
cates one of the spatial coordinates 



(W(C)) ~ e~ s (4.6) 

where W(C) is the Wilson loop around some closed contour C: W(C) = 
means that the calculation is path-ordered. S is the action 
of the worldsheet of the string. 

The LHS of (14. 6 p needs further explanation. In Euclidean field theory 
with a rectangular contour with sides of length T and L (figure 14.11) . the 
energy of a quark-antiquark pair can be found by calculating the expectation 
of the Wilson loop in the T — > oo limit [110]: 

(W(C)) ~ e~ TE ^ (4.7) 

where E(L) is the energy of the quark-antiquark pair. A good analysis of 
equation ( 14.7ft can be found in [111]. Here we summarise the result. We 
start by defining the two quark state at time t by: 

\T(t,R)) = \q(t,0)q(t,R)) (4.8) 
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Then we need to consider 



\im Q(T,R) = \im(q{T,0)q{T,R)\q{0,0)q{0,R)) 

T— >oo T— >oo 



= lim<rt(T,J2)r(0,J2)> 

1 — >oo 

n 

~ e-^ T (4.9) 



We have been able to pick out the smallest energy eigenvalue E Q in (14.91) 
because in the T — > oo limit, it will dominate all the other terms. 

All we have left to do to demonstrate (14.71) is to show that lim Q(T, R) ~ 

T — >oo 

lim W(C). However, the Wilson loop corresponds to an extended 'source' 

T — »oo 

in the fundamental representation (ie. quarks), separated by a distance R, 
being created at T = and annihilated at some general T, so we are done. 

Let us apply equation (14.61) to the Constable- Myers geometry considered 
in this chapter. The Nambu-Goto action of the fundamental string is [109]: 

S = —5— / d 2 ay/-detG M Nd a x M dbX N (4.10) 
zna J 

with Gmn denoting the string frame metric, which is related to the Einstein 
metric by a scaling factor: 

""'string c ua Einstein 

= H" 1 ' 2 f( A +Wdxl + r 2 r XI2 fi^-m^^L dw 2 

1=1 

It will be useful to recast (14.111) with an explicit 5-sphere symmetry: 



d.S 2 = H-^f^+Wdxl + i ? 2 F l/2 / (A+2^)/4!^_^ ^2 + ^ 
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We want a static solution to the equations of motion, so it is convenient to 
set cr° = t and cr 1 = x, and assume y = z = 0. We also assume that the 
string is at a constant angular position in S 5 . Then the action becomes 



S = / dt dx n(r) 



1 + m(r) 



dx 



1/2 



(4.13) 



with 



n(r) = G xx — iy _1 /2y(A+5)/4 



Equation (I4.13P is very similar to the integral found when performing Fer- 
mat's Principle for the path of light rays [102]. For large separation of 
endpoints, the trajectory which minimizes S locates itself very nearly at 
the minimum of n(r) for most of its length. Hence the complicated pro- 
cess of optimizing S is reduced to simply optimizing n(r). Furthermore, for 
small separations Ax about its minimum, S becomes 



S ~ ^^Ax (4.14) 



It is easy to show that n(r) is optimized by 



which gives 



V(Ax) = S(Ax) = Ax (4.16) 

27ra'vl — S 



with S = ^rf ■ f)4.16p shows us that this geometry gives us a linear quark- 
antiquark potential. In other words, it exhibits confinement. 
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Glueball State 

vl 1 Ltv> I ) f ill U IcllC 


Mass 

IV J.CXOO 


w 


1 nn 

l.UU 


o++* 


1.56 


Q++** 


2.08 


qH — h*** 


2.60 


qH — [-**** 


3.11 



Table 4.1: ++ glueball masses from Constable Myers geometry. Normal- 
ization is such that the ground state mass is set to one (arbitrary units). 

4.2.2 Glueballs 

The glueball spectrum for this geometry was investigated in [109]. The 
procedure for calculating holographic glueball spectra is quite straightfor- 
ward [82-84,112]: it is found by solving the equation of motion for the scalar 
field $: 

d, (v^Te- 2 ^) = (4.17) 

with the metric given in equation (14. ip . Making the ansatz $ = g(p)e lk,x 
(since the scalar field is assumed to be a small perturbation) we then look 
for solutions to (14.171) . This can be done numerically. As r — > oo, the 
metric asymptotes to AdS space, and the normalizable solution to the wave 
equation goes like r~ 4 . Hence we use the boundary conditions g(A) = 
A~ 4 ,g'(A) = — 4A -5 we then look for those solutions which obey f'{0) = 
(for continuity). We find a discrete spectrum (table |4~T|) . In other words, 
this geometry has a mass gap in the glueball sector. 

4.3 Chiral Symmetry Breaking 

The most exciting aspect of this geometry is that once we introduce quarks, 
a quark condensate forms spontaneously. This is an essential characteristic 
of QCD. 

We start by introducing quarks in the same way we did in chapter [31 
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by adding Nf probe D7 branes to the geometry. The exact dimensions 
occupied by the D7 and D3 branes are shown in table 14.21 

We make a change of variables in the two directions perpendicular to 
both branes: 

$ = r5 + ir 6 = ae w (4.18) 

and introduce one flavour of quark via a D7 brane probe in the geometry. 
The Dirac Born Infeld (DBI) action for the probe is 

Six = - {2 ^ gs J AeV-det(P[G ab ]) (4.19) 

where P indicates the pullback of the spacetime metric onto the D7 world 
volume. The pullback of a metric is defined as 



P[Gab] = %%G a , (4.20) 



with xq —>■ xg being the coordinates in the bulk, and £o ~~ * £7 being the 
coordinates on the D7 brane. We choose the static gauge, which simply 
means = & for i = 1 -> 7, and x 8 = x 8 (£o, ...,&), x 9 = x 9 (£ , ■ ■ ■ , £7). 
This action will determine how the D7 lies in the remaining t^—tq directions. 

Performing the necessary calculation we get 



(2tt) 7 « /4 ^ J s v^ 4 - fo4 / V ^ 



4.-5 

x \ll + \d P ^ + H(^±^) 8 ^W| 2 (4.21) 



with u 2 = p 2 + |$| 2 = p 2 + to 2 + u 2 . 
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X 


X 


X 


X 







Table 4.2: Choice of D3 and D7 embedding. Filled directions are marked 
with a cross. Unfilled directions by a dot. 

Applying the Euler-Lagrange equation of motion for $ we find 

(4 , 3) 

There are two important points to be noted in equation (14.221) . 

Firstly there is a an explicit $ — > e ia $ symmetry. In the massless case 
this is the U(1)a symmetry on the quarks. 



Secondly, in the UV limit (r — > 00), the equation of motion simplifies to 

*(^) =0 (4 ' 24) 

which has the simple solution $ = m + cp~ 2 . The two integration constants, 
m and c, correspond to a quark mass and a VEV for the quark bilinear qq 
respectively. 

Returning to the full equation of motion, we look for regular solutions 
that have the UV boundary condition a = m + cp~ 2 . (For the moment 9 
is turned off: it corresponds to the pion fields, which we will demonstrate 
later) . A few of the solutions are plotted in figure 14.21 

The key feature of figure l4~2l is the m = 0, c = 1.86 curve. This shows 
that this geometry demonstrates chiral symmetry breaking: there is a quark 
condensate even if the quarks are massless! This makes the Constable- Myers 
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Figure 4.2: Solutions for the ujq flow when uj^ = showing the dependence 
of the condensate on the quark mass 

geometry a very interesting model to study other properties of QCD. 

4.4 Pions and their Interactions 

We now turn our attention to the presence of Goldstone bosons. Numerical 
solutions discussed in the last section have demonstrated that this system 
dynamically breaks the classical U(1)a symmetry by forming a quark con- 
densate. We first met Goldstone's theorem [25] in chapter fl.2.81 and in that 
section we analysed the spontaneous breaking of the U(1)l x U(1)r sym- 
metry to U(l)y. We identified the 7/(958) meson as the relevant Goldstone 
boson. However, matters were complicated by the presence of instantons, 
which made the r/(958) heavier than would otherwise be expected. Here, 
at large N, there are no instantons and so the equivalent rj meson should 
still be light. 

Henceforth we will describe the pseudo- Goldstone boson resulting from 
the spontaneous breaking of the U(1)a symmetry of this geometry as a pion. 
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4.4.1 The Pion Equation Of Motion 

In the case of one flavour, the pion is simply a bound qq state, and will 
be described by the field $. The Goldstone fluctuation corresponds to os- 
cillations of the field along the vacuum manifold, which in the chiral limit 
means fluctuations in the angular 9 direction of $. To leading order, if we 
have a background configuration for the D7 brane described by o"o, we can 
look at small fluctuations in the 9 direction. 



The action for these fluctuations (to quadratic order) is given by ex- 
panding the DBI action (14. 191) : 



2tt 2 



(2vr) 7 «' 4 ^ 



dp d A x e^R A Q^\ + (d p a o y 



i + 1 g"W W + l 9^9eed u 9d^ 



2 1 + (d p a o y 2 1 + (<V 



And the resulting Euler-Lagrange equation of motion is: 



MW£ tr ^ (p 2 + a 2 ) 2 + ^ (1 " 3)/2 (p 2 + ^) 2 -l 



(7, 



d ( e *g 2 \ = 

dp\^l + {d p alY ° pJ ) 

We look for values of MR as a function of m that give regular solutions to 
(14.26p . The UV boundary condition we impose is / = 1/p 2 , which reflects 
the fact that pion has the same UV scaling dimension as qq. The results are 
plotted in figure I4T31 There is a massless pion when m q = 0, in accordance 
with Goldstone's theorem. 
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Figure 4.3: Quark mass versus meson mass, showing the massless pion when 
m q = 0. We plot the square root of the quark mass to demonstrate the linear 
relationship. 



4.4.2 The Vector Meson Equation Of Motion 

There is an additional field present on the D7 brane which we have left 
undiscussed: it is the gauge field partner of the scalar $. The field strength 
of this gauge field enters in the standard way as 2Tta'F ab in the square root. 
The leading Lagrangian for this field, in a background configuration o"o is 



fill 



There is an additional term that could be added to the DBI action: the 
Wess-Zumino term, which gives the coupling of the four-form to the 
gauge fields. It is not included here because when we calculate the equations 
of motion for the gauge fields, this term is only relevant for the gauge fields 
with a vector index on the S 3 . We are only interested in states that carry 
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n 



AdS case 



CM case 




1 

2 
3 
4 



2.83 
4.90 
6.93 
8.94 
11.0 



2.16 
4.85 
7.05 
9.20 
11.3 



Table 4.3: Vector meson spectrum comparing CM and pure AdS back- 
grounds 

no SO (4) R-charge. We write the gauge field as 



There are two asymptotic solutions to equation (14.291) . The simplest is 
M — 0, g(p) = constant. This corresponds to introducing a background 
gauge field associated with U(l) baryon number in the field theory. It is of 
little interest to us here. 

The second solution is g(p) ~ 1/p 2 , and this has the right dimension 
and symmetries to be identified as dual to the operator q / ~f fl q. By seeking 
smooth solutions to the equation of motion we can determine the vector 
mass spectrum. The results are shown in table 14.31 and compared to the 
equivalent result in the pure AdS$ x S 5 spacetime, which we calculated in 
section 13.2.21 



= g(p) sm(kx)e fl 



(4.28) 



The resulting Euler-Lagrange equation of motion is then 



+ d p 




(4.29) 
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• SU(2) 
■ SU(3) 

♦ SU(4) 
i SU(6) 
— N^inf. 



0.1 0.2 0.3 0.4 0.5 0.6 

Figure 4.4: A plot of m p vs in the Constable-Myers background on the 
left. Lattice data [1] (preliminary, quenched and at finite spacing) for the 
same quantity is also shown on the right. 

4.5 m p vs. ml 

We can compare the dependence of the rho meson mass on the pion mass 
squared. This is an interesting thing to do because the same result has been 
computed for large N in lattice simulations [1], so a direct comparison of 
gauge/gravity and lattice results is possible. We show the result in figure 

We have chosen units in the Const able- Myers geometry such that we 
can compare directly with the lattice results of [1]. It is striking that not 
only does the lattice data display the same linearity as the gauge/gravity 
result, but that the slopes are so similar (0.52 and 0.57 respectively). This 
is one of many examples where AdS predictions match QCD much better 
than one would naively expect. In the next section we explore some more 
predictions. 
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Chapter 5 
AdS-QCD 



5.1 Premise 

In 2005 two very interesting papers [2,113] were published that for the first 
time dispensed with trying to deform supergravity solutions to mimic QCD. 
Instead the authors aim was to "start from QCD and attempt to construct 
its five-dimensional holographic dual". Since QCD has an infinite number 
of operators, they had to be careful about which operators they were going 
to try and include. [2] focused on holographically modelling the dynamics 
of chiral symmetry breaking in QCD. To this end, we would expect only a 
small number of operators to be influential. These would be 

• the left-handed and right-handed currents corresponding to the SU(Nf)i x 
SU(Nf)n chiral flavour symmetry 

• the chiral order parameter 

So in keeping with the unique operator-field matching, this would require 
three fields on the five-dimensional holographic side. Their details are given 
in table 15.11 
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AD:0(x) 5D:0(x,r) p A mjj 

q L rt a q L A%l 1~~3 0~ 

q Rl H a q R A Rlx 13 

q a R ql (2/z)X<*P 3 -3 



Table 5.1: Operators/fields of the model of [2] 



The choice of for each field is determined by the relation (A— p)(A + 
p — 4) = m 2 which we discussed in section 12.71 The metric, also an ansatz, 
is simply AdS$, which is the simplest choice we could make. In keeping 
with [2], in this section we will write the fifth dimension in terms of z, 
instead of the usual r. The two are related by z = r -1 . Hence z ~ is 
the ultraviolet part of the spacetime, and z ~ z m is the infrared of the 
spacetime. 

ds 2 = -z~ 2 dz 2 + z~ 2 rj^dx^dx u (5.1) 

The use of an AdS geometry implies that the dual theory will be conformal 
- not something we want in a dual to QCD. So to give the theory a mass 
gap, the spacetime is defined to end at some IR scale z = z m . In the window 
< z ^ z m , the background gauge is conformal, and the coupling doesn't 
run. 

The 5D action is 

S = Jd 5 x^Tr{\DX\ 2 + 3\X\ 2 --^(F 2 + F 2 

where D^X = d^X - iA Lfl X + iXA R ^ A LiR = A a L R t a , and = d„A„ - 
d u A fM — i[A^, A U ].X is a Higgs-like field, holographically dual to the quark 
mass and condensate. Its value in the UV will be: 

X (e) — ► 2m q e + 2e 3 a as e — > (5.3) 
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where m q is the quark mass matrix, and a is the quark condensate matrix. 
Both are diagonal matrices. 

There are now four free parameters: z m ,m q ,a, g§. g§ can be fixed by 
operator product expansion for the product of currents, and we come to this 
a little later. We also take this opportunity to focus on just two flavours of 
quark, so in table I5TT1 a, f3 = 1,2; a,b = 1,2,3; and t a = <J a /2, where a a 
are the Pauli matrices. This model with three flavours of quark have been 
investigated [114], and the results are good. Here though, the third flavour 
would provide an unnecessary complication. Models beyond three flavours 
would be futile because the lightest three quark flavours are of roughly the 
same mass and thus have an approximate SU(3) symmetry. The charm 
quark however is much heavier than the previous three quarks, and so the 
SU(4) symmetry is badly broken. 

As with the AdS / CFT correspondence, we claim that precise correspon- 
dence is obtained by equating the generating functional of the connected 
correlators in the four-dimensional theory with the effective action of the 
five-dimensional theory, with the UV boundary values of the 5D bulk fields 
set to the value of the sources in the 4D theory: 




(5.4) 
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5.2 Fixing the g$ coupling 

The idea in this section is to calculate the vector current two-point function 
in the holographic scheme, and then by equating the result to that obtained 
from perturbative QCD we can hence fix g§. 

We introduce the vector field as V = (A L + A R )/2, choose the gauge 
V r (x, r) = 0, and by using the Euler-Lagrange equation on the holographic 
action (15.21) . we obtain an equation of motion: 



[dJ±d z V2(q,z)^+^V«(q,z)] x = 



(5.5) 



We then substitute (I5.5P back into the action (15.21) . and we are left with the 
boundary term 

1 r . /I \ 

(5.6) 



S = - A / d A x [ -V?d 7 V» a 



2g!J \z " 

So, by the claim made above (15 .4p . equation (15.61) is the generating func- 
tional for QCD. To obtain the vector current two-point function, we simply 
functionally differentiate (15.61) twice with respect to the source Vo, with Vo 
defined in V^(q, z) = V(q, z)VQ(q). iy^^q) is the Fourier transform of the 
source of the vector current J® = q^^q). Doing the necessary maths, we 
get 

j d*x(J«(x)J b u (0)} = 5 ab (q,q u -q 2 g, u )ll v (Q 2 ) (5.7) 
tt / 2\ 1 d z V(q, e) 

IM-? ) = (5.8) 



For large Euclidean Q 2 (= —q 2 ) we only need to know V(q,r) near the 
boundary, 



V{Q, z) = l + ^f- ln(Q 2 z 2 ) + ... (5.9) 
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which upto contact terms gives 



iiv{Q 2 ) = -^- 2 HQ 2 ) (s.io) 



We can calculate exactly the same quantity using Feynman diagrams [115, 
116]. The leading-order diagram is the quark bubble 



ny(Q 2 ) = -^ln(Q 2 ) (5.11) 



And so comparing (15. 10p and ( 15. Ill) we can fix g§ as 



, 12tt 2 

= ivr (5 ' 12) 



A cautionary note 

The matching we have just performed is a little naive. One should match 
the gravitational theory to QCD only at the point where the QCD coupling 
becomes non-perturbative, for this is the point where the gravitational dual 
will be weakly coupled. And of course, when the QCD coupling is non- 
perturbative, gluonic contributions become important. In all of the analysis 
to follow, we keep g$ fixed as in equation ( 15.121) . However we have performed 
an analysis that treats g$ as a free parameter: this indicated that the optimal 
value of #5 is 5.19, rather than the 3.54 that (15.12ft suggests: this is a 
discrepancy of 32%! So for the rest of this section we add the disclaimer 
that non-perturbative effects could have a significant effect on our results. 
We also remind the reader that AdS/CFT correspondence is a large N c , 
large 't Hooft limit (gy M N e A > 1). In calculating g 5 we have just 
assumed N c = 3; this is a further source of error. We return to the matter 
of g$ in chapter [6] 
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5.3 Calculations 



5.3.1 Vector Mesons 

There are now only three parameters in our model: z m ,m q and a. The 
hadrons of QCD correspond to the normalizable modes of the 5D fields. 
Requiring the modes to be normalizable means that when we substitute the 
chosen modes back into the 5D action, the action remains finite. As usual, 
this means that the solution must die away sufficiently rapidly in the UV, 
and be smooth over the IR boundary (here given by z = z m ). This gives us 
two explicit conditions for a general solution ip(z). 

ip(0) = d z Tjj(z m ) = (5.13) 

There is also a third constraint: that of orthogonality. ip n and ip m must be 
orthogonal to one another if they are to represent distinct mesons. So we 
also have f [dz/ z)^ 2 = 1. Two constraints is sufficient to solve a second or- 
der differential equation: having three constraints guarantees us a discrete 
spectrum. This is clearly a pre-requisite for a sensible model of QCD. 



A simple check shows that 



, ) = £jM£)^£l (5.14) 
' q z — mi + te 

o r 



is the Green's function to ( 15. 51) . It can also be shown that V(q, z') of equa- 
tion (15.81) is given by — (l/z)d z G(q; z, z') at z = e. Putting this together 
with equation ( 15.81) we find 

I fh> 4 E - (5.15) 
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which on comparison with equation II. 381 of section [1 .2.81 in the limit e — ► 0, 
allows us to extract the decay constant F p : 



F 2 P = -Ami*Y 

95 



(5.16) 



5.3.2 Axial Vector Mesons 

Our analysis of the axial vector mesons closely follows that of the vector 
mesons. However things are complicated by the mixing of the pion field ir 
and the longitudinal axial gauge field <j). We start with the action, expanded 
to quadratic order: 



S 



v(zf 



2z A 



(5.17) 



where we have defined v(z) = m q z + az 2 ,A = (Al — Ar)/2 and X = 
X exp(i2n a t a ). We make the gauge choice A z = and A^ = A^± + <9 M 0. 
Then the equations of motion become 



( -cU« ) + Q -j 



9W . a 



— 5-(7T - 

z 6 

9 2 5 v 2 



d z ir a 









(5.18) 
(5.19) 
(5.20) 



The di meson is a spin-1 particle, and so is the solution to equation (15.181) 
with the same boundary conditions we imposed for the vector mesons, 
namely 



<Mo) = o 



d z ip ai {z r , 







(5.21) 
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The normalization condition is also the same: J(dz/z)xjj 2 1 = 1. Similarly, 
the decay constant F ai is given by (I5.16p . but with p replaced by a\. 



< = H? «(<0A) 2 (5-22) 
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The pion can be found by solving (I5.19P and (I5.20p simultaneously, subject 
to the boundary conditions <f)'(z m ) = 0,0(0) = and 7r(0) = 0. The pion 
decay constant can be found by considering the axial current correlator of 
(JESS): 

(01^(0)1^))= z/^ (5.23) 

According to (15.231) . n^( — q 2 ) will have a pole at q 2 = m 2 . In the limit m n = 
0, rU(— g 2 ) ~ —f 2 /q 2 . Hence, on comparison with ( 15.1 6ft . we conclude 



SI - A 9 -^ (5-24) 

95 e 



5.4 Results 

From this simple model we can predict six quantitiej*]: the mass and decay 
constants of the p mesons, the pions, and the a± mesons. We still have 
three inputs: the quark condensate, the quark mass, and the position of 
the infrared cutoff. Performing a best-fit analysis gives a rms error of 9 
%.(Note e rms = ^2 ((50 / 'O) 2 / 'n) 1 / 2 with O the observable and n equal to 
the number of predictions minus the number of inputs, so here n = 3). The 

results are displayed in table 15.21 

*in [2] they claim to predict seven quantities; the seventh being the n — p coupling. 
(7 p7r7r . However this term would presumably receive an important contribution from F 3 
terms in the holographic action, which haven't been included. It could be argued, though, 
that a good agreement between the prediction for g p7r7r and its experimental value in the 
absence of F 3 terms indicates that these terms have little influence 
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V. ./ y_r Ovl V < 4 YJ L V 


Measured fMeV") 

i-V J-V^CliO LAX CU I IVlv V / 


Best Fit fMeV") 


777, _ 


140 


141 


m„ 


776 


832 


m ai 


1230 


1220 


u 


92 


84 




345 


353 


Fii 2 


433 


440 



Table 5.2: Results of model from [2] for QCD observables. m q = 2.3MeV, 
a = (327MeV) 3 , z m = 323(MeV)~ 1 

5.5 Improvements 

Since the publication of [2] in December 2005, there have been many further 
publications, suggesting improvements to the simple model presented in the 
preceding parts of this chapter. We consider here one of the more well- 
known extensions, achieved by including a dilaton in the action (15.21) and 
generalizing the metric (15.11) to: 

S = Jd 5 xe-*W v^Tr||M| 2 + 3|X| 2 -^2(F| + F|)J(5.25) 
ds 2 = e 2A{z) (dz 2 + VlMU dx»dx") (5.26) 

As yet, $>(z) and A(z) are undetermined. We can constrain them with the 
following pieces of physics: 

1. The geometry must return to AdS in the UV (z — > e), so the linear 
combination $ — A must asymptote to log z for small z 

2. A simple flux model [117] shows that meson masses obey the relation 
rn? n ~ n for each radial excitation n. This agrees with experiment 
(figure I5~TT) . To reproduce this phenomenon, $ — A must asymptote 
to z 2 for large z. It is worth noting that the model [2] we analysed in 
the previous section predicts a meson mass spectrum with m 2 ~ n 2 
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Figure 5.1: Experimentally, the p mesons follow very closely the relationship 
w? n ~ n. The dots are the physical values 

There are still an infinite number of choices we could make for A(z) and 
$(z) that would obey these constraints. We make the simplest choice, and 
set A(z) = — \og(z) and $ = z 2 . Hence the action and metric become 



Now we repeat the procedure of section 15.31 to find the meson masses and 
decay constants. 

5.5.1 Vector Mesons 

Repeating the process of section 15.3.11 the equation of motion becomes 



As usual, we look for normalizable solutions to ( 15. 29ft . However, now our 
choice of boundary conditions is much simpler than when there was no 
dilaton term. There is no hard IR cutoff at there was previously) , 

and the action is defined for all z. So the ambiguity of (15. 131) is removed 
when we have a theory with a smooth wall cutoff (although we have re- 



ds 2 



S 





(5.29) 
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introduced some ambiguity through our choice of A(z) and <&(z)). Requiring 
that the solution is simply normalizable gives us a discrete spectrum: 



v n (z) = z\\ ^-rLKz 2 ) (5.30) 
n + 1 



where L™ are the associated Laguerre polynomials. The masses and decay 
constants are found using the same formulae in section l5.3.1l the only differ- 
ence being the presence of the non-zero dilaton, which changes the results 
to: 

m 2 n = 4(n + l) (5.31) 
Fl = ^ (5.32) 

Note that the mass spectrum goes like m\ ~ n, as we set out to achieve. 
5.5.2 Axial Vector Mesons 

The treatment of the axial vector mesons follows that of section 15.3.21 We 
define the axial gauge field as A = | (A L — A R ), and the equation of motion 
is found to be 



d z [ ^d z a n J + (m 2 n - ^X(z)A ^a n = (5.33) 



with the linearized equation of motion for the Higgs-like field X(z) given 
by 

' -z 2 \ -z 2 

d z ( e —d z X(z) + 3^-X(z) = (5.34) 

Z J Z 
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and, of course, in the UV the Higgs-like field must asymptote to the usual 
holographic quark mass (m q ) and condensate (a) terms: 

X(z)^~m q z+±az 3 (5.35) 

Unfortunately, although it is not yet clear, including the dilaton term in the 
action (15.271) ruins the predictability of the axial vector meson sector. We 
demonstrate why in the next section. 

Loss of predictability 

For large z, the equation of motion for X(z) (15.341) becomes 

X" - 2zX' + —X = (z > 1) (5.36) 

(15. 36ft has two linearly independent solutions. One asymptotes to e z in the 
IR, whilst the other tends to e~°' 75z 2 . The solution must be finite for all 
z, so we must discard the first solution, and hence conclude that m q and 
a are proportional to one another. This is not what one wants in a theory 
with anomalous symmetry breaking. We conclude that including a dilaton 
term $(z) = —z 2 greatly improves the vector meson sector, but removes 
anomalous symmetry breaking from the holographic field theory. 

Further extensions 

We only included terms quadratic in X in the holographic action ( 15. 27ft . 
The previous section showed that this is not sufficient to model the axial 
vector meson well. Including such terms would introduce nonlinearity into 
equation ( 15.341) . This in turn would make the relationship between m q and a 
nonlinear: an essential result if we are to holographically model spontaneous 
symmetry breaking. This avenue is left to other authors. 
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Many other aspects of QCD phenomenology have been successfully broached 
by AdS-QCD models. Baryons, strange quarks and higher spin mesons have 
all been addressed [114,118,119]. Four-point current-current correlators rel- 
evant to the AI = 1/2 rule and the Bk parameter for K- meson mixing are 
analyzed in [120,121]. Heavy quark potentials are computed in [122]. The 
AdS/QCD model is related to light cone QCD in [123,124] allowing form 
factor computations. Form factors for mesons can also be found in [125,126]. 

Properties of QCD at high temperature and density and the deconfme- 
ment transition have been analyzed in this context in [127-133]. 

Such models have also been adapted to describe walking [134] techni- 
colour [135,136] dynamics for electroweak symmetry breaking in [137-141]. 



84 



Chapter 6 

Improving the IR of 
Holographic Descriptions of 
QCD 

As we saw chapter^ a surprisingly good description of QCD can be obtained 
from a simple five dimensional gauge theory on a truncated AdS space [2]. 
However, there are many short-comings of the model: 

• The use of an AdS geometry implicitly means that the background 
gauge configuration is conformal 

• Hence the UV does not become asymptotically free 

• The existence of a mass gap is imposed by hand by the introduction 
of an infrared cutoff. It is not the result of a running coupling 

• The fields that holographically describe the quark bilinears are in- 
cluded phenomenologically 

• The solution for the field holographic to the quark mass and conden- 
sate is put in by hand - it is not dynamically determined by either the 
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gauge configuration or the quark mass. 

• The excited meson mass spectrum scales like the excitation number 
n. Experimentally, the masses scale like y/n, a result which a simple 
flux model can reproduce (see [117]) 

It would be useful to address the above points, and some progress has 
been made [117, 142, 143]. Here we take the opportunity to examine the ap- 
proach taken in [142]. The idea was to take the well-documented Constable- 
Myers supergravity solution of chapter HJ whose dual describes many of 
the characteristics of QCD, most crucially, chiral symmetry breaking. The 
salient features of the model are 

• The background gauge configuration in which the quarks live is non- 
supersymmetric, and has a running coupling 

• The mass gap arises naturally as a result of the non-supersymmetric 
gauge configuration 

• The holographic dual of the quark bilinear is explicit in the string 
configuration 

• The quark condensate is a prediction of the gauge configuration, and 
is determined by the quark mass 

These points go a considerable way towards addressing the short-comings 
mentioned earlier. We will, however, continue to adopt the phenomenolog- 
ical approach of treating the background as describing an N = 3 theory 
rather than N = oo. Furthermore, the string theory construction can only 
realise a U(l) axial symmetry and does not provide a holographic descrip- 
tion of the axial vector mesons. We artificially include the appropriate fields 
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to provide a non-Abelian chiral symmetry and axial vector states. This is 
in keeping with the phenomeno logical spirit of [2, 113]. 

In this chapter we compute with our phenomenological model the masses 
and decay constants for the pion and the rho and a\ vector mesons. We 
find that the model gives comparable predictions to the pure AdS models 
within 12% of the QCD values. We believe these results provide support 
for the robustness of the predictions of these holographic models. 

The geometry we propose returns to pure AdS space in the ultra-violet, 
so we do not address here the absence of asymptotic freedom in the grav- 
ity description. As pointed out in [144], the gravity theory should only be 
used up to a UV cut off, corresponding to the scale at which QCD switches 
from perturbative to non-perturbative behaviour. Above that cut off the 
gravitational dynamics must become non-perturbative with its loop correc- 
tions completely dominating the classical results. The correct UV dynamics 
should be encoded at that cut off by correcting the values of higher dimen- 
sion operator couplings. In principle, these can be tuned in the AdS/CFT 
approach to produce the holographic equivalent of a perfect lattice action. 

As a small example of these ideas we consider the matching of the five 
dimensional gauge coupling in the UV. In [2, 113] this coupling is matched 
to the perturbative result for the vector vector correlator in QCD. The 
AdS gravitational dual presumably describes a strongly coupled conformal 
theory in the UV and so the correlator behaviour matches the logarithmic 
result of the conformal but weakly coupled UV behaviour of QCD. It is 
surprising that the numerical coefficient of the log term can be matched 
though. In chapter we indicated that this matching might underpredict 
#5 by as much as 30%. We perform a similar analysis here, and find that 
once again the ideal value of g§ is higher than that indicated in [2, 113], 
this time by 18%. This provides a measure of non-perturbative corrections 
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at the scale of matching to the strongly coupled regime of QCD. We leave 
attempts to further improve the UV of the theory to chapter [7J 

Finally, as we showed in section 15.51 an appropriate change to the IR 
behaviour of the dilaton can correct the n scaling of the tower of excited 
p meson states. We have tested our model in this respect but find only a 
marginal improvement over the pure AdS case. This is a sign that, although 
our geometry describes a non-supersymmetric gauge configuration, it is still 
not a perfect description of QCD and work remains to be done on improving 
the geometric background. 

6.1 The generic 5D holographic model 

For a generic five dimensional holographic model of QCD, we can write the 
action as 



with DfJJ = dJJ - iA L pU + iUA Rfl . The field U(x,r) = exp(in a (x, r)T a ) 
describes the pions produced by the breaking of a U(Nf) chiral symmetry 
with generators T a . The non-Abelian gauge fields Al and Ar couple by 
left and right action on U and holographically describe the axial and vector 
mesons (we describe this in more detail later). The a is a Higgs-like field and 
we only consider fluctuations in its r direction. It holographically describes 
the quark mass and condensate. A non-zero value for this field will break 
the U{Nj)l x U(Nf)n symmetry of the action down to U(Nf)y- 

If we set 4> = and C a = 3|X| 2 , we reproduce the model discussed in 
chapter [5j 

However, in this chapter, we are going to use a compactified model of 



S 




(6.1) 



SS 



the Constable- Myers solution, presented in chapter HI 



6.1.1 Truncating the Constable-Myers Geometry 

The truncation we propose is incredibly simple. For convenience we repro- 
duce the equations of motion (14.261) . (14.291) for the eight dimensional fields 
defined on the surface of the D7 brane. The functions /„ represents the 
pion fields, and g n the vector fields. 



R 2 Mlet>g ( { P 2 + al? + l \ {1 - 5)/2 (P 2 + alf - 1 

>/i + ( w V(p 2 +^o 2 ) 2 -i; (p 2 +^o 2 ) 2 ° 0Jn 



2 



e*g^/l + (d p cr yM*R 2 g n (p)H 

+ d P ( h = =dMp) /r^^ TTT ) = (6 - 3) 

Q is defined as in ( 14.231) . The reader is reminded that the metric is also 
defined in eight dimensions: 

ub 8D — 

(6.4) 

Since we want to propose a theory that is dual to QCD, the fields must 
not have any components on the three sphere. The simplest way to do this 
is to claim that all of the fields f n , g n , <fi and a are constants on the 3-sphere. 
They are only functions of r and x only. 



W + T) 
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On making this choice, we can then either work directly with (16.21) and 
(16.31) . whilst remembering that 0, n and a are constants under the SO(3) 
isometry, or we can explicitly recast the entire ensemble in five dimensions. 
This is how it was presented in [142], and we reproduce the five dimensional 
theory here. 

The five dimensional metric is 

ds 2 = H^ 2 fV 8 V^dx^dx" + H^f^hdr 2 (6.5) 

i=0 

with 



(a(r) 2 + r 2 ) 2 + 1 (a(r) 2 + r 2 ) 2 - 1 1/2 _ 

7 ((7(r) 2 + r 2 ) 2 -!' (a(r) 2 + r 2 ) 2 ' 7 



(6.6) 



The explicit values of 0, g 5 and £ CT are 



e = if 5/4 / 15/16-V39/4^/2 r 3 (1+d 2 ) -l/2 ( g y) 

g\ = 4^2^1/2^/8+^/4^ (6g) 

^ = ^f V ^ 9/4 9TVT+^ (6.9) 

So on substituting (16.91) into (16. ip we get the five dimensional action. Note 
that as r — > oo, the entire model reverts back to the simple AdS model 
considered in section [51 

We now go on to calculate the masses and decay constant of the vector 
and axial meson sectors. 
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6.2 Vector and Axial Mesons 



6.2.1 Vector sector 

In exactly the same way that we calculated the vector and axial meson 
masses and decay constants in chapter [51 we repeat here. 



Just as in section 15.3.11 we look for solutions to the vector equation 
of motion that are of the form V°(x,r) = V£(r) exp(iqx) with the gauge 
V r a {x,r) = 0. We find 



with K x = f^hr^l + a 2 )- 1 ' 2 and K 2 = if/ 1_(5 / 2 /i 2 r 3 (l + a 2 )- 1 ' 2 



The p mesons are interpreted as the normalizable modes of (I6.10p . mean- 
ing that we impose the boundary conditions V'(0) = to ensure the smooth- 
ness of the solution, and V(oo) = to ensure the normalization is finite. 

Our analysis for the decay constants follows the same path as that laid 
out in section l5.3.1l For large N, the vector current correlator can be written 
as the sum over p resonances: 



In order to find F p , we proceed by finding the Green's function solution to 



d^K^drV^r)) + q 2 K 2 {r)V«{r) = 



(6.10) 




(6.11) 



(16.1 Op . Imposing the completeness relation K 2 {r)%lj p {r)ip p {r') = 5(r — r') 
on the set of eigenf unctions one finds 




(6.12) 
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The generalisation of equation (15.81) is 



Ki(r)d r v(q, r) 



gh 2 



(6.13) 



Using the fact v(q,r') = [Ki(r)d r G(q; r, r')] r=00 , we find 



1 (KArWJr)) 2 
Kv(-q 2 ) = ~ Hm ~ 2 V \ £ \ (6.14) 
™ 91 ^ {q 2 - m 2 p )m 2 p 

Comparing this to (16.111) we identify the rho decay constant as: 



F 2 p = limi^M^r)) 2 (6.15) 



6.2.2 Axial sector 

The axial vector field is defined as A a p = (A a Lp — A a Rp ). We choose the gauge 
such that A ar (x,r) = and A% = A^ ± + <9 M 0. There are three equations of 
motion: one for each of the a\ meson field, the 7r field and the longitudinal 
axial gauge field <p: 

[d r (K x (r)d r Al(r)) + q 2 K 2 (r)A;(r) - g 2 5 a (rf K 3 (r) A%rj\ ± = 0(6.16) 
9 r (^ 1 (r)9 r 0)+^ 5 V(r) 2 ^ 3 (r)(7r a -0 a ) = 0(6.17) 
-q 2 K 1 (r)d r (f) + g 5 2 K 4 (r)a(r) 2 d r 7r = 0(6.18) 

where K 3 (r) = Hf 3 / 2 - s / 2+A / 2 h 3 r 3 (l + a 2 )- 1 / 2 and K 4 (r) = f 1+A / 2 h 2 r 3 {l + 
d 2 ) -1 / 2 , and we have separated A± as A i± (x,r) = A(q,r) exp(iqx). 

Looking for the normalizable solutions to (16.161) requires us to set the 
boundary conditions ^ a i(oo) = and d r ip a i(0) = 0. We can then predict 
the a>i meson mass spectrum. And in an almost identical analysis to that 
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V/ UOVl V < \ V J L\. 


Measured fMeV") 


New Model fMe\H 


Previous Model fMeVl 


Tfl-rr 


140 


139 


141 


m 


776 


743 
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m ai 


1230 


1337 


1220 


u 


92 


84 


84 


b P 


345 


297 


353 


-Tat 


433 


491 


440 



Table 6.1: Comparing the best fit results in [2] and the model presented 
here 

we performed for the vector sector, the a\ decay constants are given by 



F 2 ai = lim-^W^M) 2 (6.19) 

Looking at equations (16.1 7p and (I6.18P we see that the pion and longitudinal 
gauge fields mix. This makes finding numerical solutions to (16. 17ft and (16.181) 
very hard. In principle we'd have to search for solutions with two unknowns: 
the pion mass, — q 2 , and the ratio of the <fi and ir fields at r = 0. So instead, 
we use the Gell-Mann-Oakes-Renner relation {m^f 2 = 2m q c) [145], to fix 
m n . Here we introduce this relation by hand, but it is possible to show 
that in the geometry studied in chapter [51 the equation can be derived 
explicitly [2]. This relation eliminates one of our free parameters, and makes 
the numerical search much easier. 

This model has two free parameters: m and b. m represents the quark 
mass (this is a model with only two flavours of quark, which we have as- 
sumed to be of equal mass), and b, which in the original supergravity solu- 
tion represented the size of the conformal symmetry breaking. Here then it 
roughly corresponds to Aqcd- Note that this model has the same number 
of free parameters as real QCD. 
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6.2.3 Results 

The results of the model are displayed in Table 16.11 We compute six QCD 
meson parameters for our fits (we do not include g pn7T ). Our model has two 
free parameters (after fixing g§ phenomenologically as discussed above), b 
corresponding roughly to the strong coupling scale A, and m corresponding 
to the light quark mass. The model therefore has the same number of free 
parameters as real QCD. 

In the first model, A, we match b and m by demanding that we correctly 
reproduce m n and m p . In order to do this, we must set At = 264.5 MeV 
and m = 2.16 MeV. This gives a prediction of 325.8 MeV for the scale of 
the quark condensate. The overall rms error for this model is 12.8%. For 
comparison we also reproduce the pure AdS fit to the same parameters found 
in [2]. That model has three free parameters, the value of the IR cut off, 
the quark mass and the quark condensate and is therefore less predictive. 

In model B, we perform a global fit to all observables. This gives A& = 
253.2 MeV and m = 2.24 MeV, with the characteristic scale for the quark 
condensate 311.9 MeV. The overall rms error for this model is 11.6%. Again 
we reproduce the equivalent pure AdS model fit for comparison. 

It is also possible to calculate the rho-pion-pion coupling constant in 
this geometry [142]. This was also done in [2]. However, as we noted in 
section 15.41 this term would receive important contributions from F 3 terms 
in the holographic action, which we have thus far neglected. Nevertheless, 
for comparison with [2] we do perform this result, and find that in this 
geometry g pn7T = 4.81 MeV. This compares well with the experimental value 
of 6.03 ±0.07 MeV [24]. 
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g 5 coupling 

In section 15.21 we pointed out that the perturbative matching of g$ may not 
be ideal. So once again, we have performed a fit with g 5 being considered as 
a extra free parameter. This indicated that the optimal value of g§ is 4.36, 
which is 18% higher than 3.54, the value suggested by the perturbative- 
matching result in equation (15.121) . 

6.2.4 Conclusions 

We have adapted a string theoretic model of chiral symmetry breaking to 
a phenomenological description of QCD. The model we have proposed goes 
some way towards addressing the inconsistencies of simple AdS slice holo- 
graphic QCD models [2,113]. The background geometry of our model is non- 
supersymmetric, and it is the smooth variation of this geometry with the 
radial direction r that provides a mass gap, without the need for an artificial 
hard IR cut-off. In addition, the dual field to the quark mass/condensate 
operator is a natural part of the geometrical set-up with the value of the 
condensate being determined by the quark mass. 

However, this is still a phenomenological approach in that we introduce 
extra fields and symmetries by hand into the model in order to describe 
the full pion and axial vector sectors. Formally there is no geometric string 
interpretation for this system. We also treat the background as though it 
describes an N = 3 rather than an N = oo field theory by matching the 5D 
gauge coupling to QCD. 

We find that the predictions of this model match experimental results to 
within 12%. This model is a little more predictive than the pure AdS slice 
models since the condensate is dynamically determined by the geometry. 
The best fit is in fact a few percent worse than the AdS slice models but 
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hopefully the theoretical improvements represent at least a moral victory. 
In any case one would naively have expected errors of order a few 100% in 
all of these models so the closeness to QCD across a range of holographic 
models supports the robustness of the approach. 

A drawback of these models to date has been that the geometry returns 
to AdS for large r, meaning that the field theory is not asymptotically 
free in the UV. Incorrect physics in the UV may affect the strong coupling 
regime in the IR [144]. Here we investigated corrections to the matching 
of the 5D gauge coupling to naive perturbative QCD results. We found 
that this coupling's value should be changed at the 18% level indicating the 
size of non-perturbative effects. In the future one might hope to study the 
importance of higher dimension operators in the IR physics as well. 
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Chapter 7 

Improving the UV of 
Holographic Descriptions of 
QCD 

This thesis has investigated the possibility of finding a holographic descrip- 
tion of QCD. However, QCD is asymptotically free and thus its dual is 
strongly coupled in the UV. We don't know how to solve string theory in 
the strongly coupled regime, and so a naive step is to simply eliminate this 
part of the duality by cutting off the gravity dual in the UV at some as yet 
undetermined energy A. The dual field theory is then only defined upto A. 

In this chapter we investigate the effect of including such a sharp UV 
cutoff in the gravity dual. The first task is to identify operator-field matches 
at r = A rather than at r = oo, which we discuss next. 
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7.1 Operator-field matching 

The matching of operator dimensions to field solutions is complicated by 
the presence of a UV cutoff, A. Without it the matching is performed at 
r = oo: now the matching must be performed at r = A. This makes things 
more difficult for two reasons: 

1. As we discussed in section II. 2. 6[ fields can have anomalous dimensions 
in the bulk which become small as r — > oo, and so don't require 
consideration when matching at r = oo. By introducing a finite cutoff 
r = A, we need to match not to the classical dimension, but to the 
anomalously- adjusted dimension of the fields. 

2. The presence of large couplings for higher dimension operators before 
QCD can be matched to the perturbative gravity theory would be 
another signal of non-perturbative phenomena. Formally there are an 
infinite set of such couplings. They are called irrelevant perturbations 
(because they become small in the perturbative regime), and will ap- 
pear in the gravitational dual as deformations of the metric which 
grow at large r. 

These two points seem almost unsurmountable, but all is not lost. There 
are an infinite number of operators in QCD anyway, yet we have shown in 
chapters [MS that including just a handful gives us very close matches to 
experimental values. In addition, we do have some idea of what to expect 
from lattice QCD results from a subject called 'perfect actions' [146]. 

For the rest of this section, we first see the effect of imposing a hard 
UV cutoff, without worrying about the effects of anomalous dimensions or 
higher dimension operators. The results are very good. We then consider 
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the possible effect of anomalous dimensions and higher dimension operators, 
and find that for the cases we consider, the results are changed minimally. 



7.2 A UV Cutoff 



Here our starting point is the model described in section 15.51 chosen both 
because of its simplicity, and because the model is smooth for the entire IR 
spectrum. Whereas before we were looking for solutions to the equations 
of motion that went like f n { r ) ~ r ~ 2 as r — > oo, we now look for solutions 
that go like f n {r) ~ at r = A. A is the cutoff, above which QCD be- 
comes perturbative or, in the holographic dual, the scale below which the 
supergravity approximation is trustable. A will be determined phenomeno- 
logically. For convenience, we reproduce the equation of motion for the 
vector meson spectrum (15.291) here: 



The results for the first five radial excited states with several different values 
of A are shown in figure I7TT1 Note that the A = oo reproduces that of section 
15.51 as expected. 

The best phenomenological fit is when the cutoff is set at 194 MeV, and 
gives an rms error of just 1.8%. This is a startling result, not just because 
the fit is so good, but also because it suggests that the regime of validity 
of a gravity dual to QCD is rather small. Similar conclusions were made 
in [144]. 

It is worth pointing out that with this value of A, we predict the next 3 
excited p mesons to be at 2320 MeV, 2475 MeV and 2626 MeV. Experimen- 
tal searches have reached up to 2510 MeV, and so far the highest excited p 
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Figure 7.1: The mass squared of the n lightest p meson excitations for 
different values of the UV cutoff. The dots represent the experimental 
values. 

meson found is the final one listed above at 2150 MeV. 

7.3 Decay constants 

The decay constants F p can be found by substituting the regular solutions 
fn{r)p^{x) back into the five dimensional action, and integrating over r, in 
exactly the same manner as section 15.3. 1L Repeating the result here for 
convenience, the decay constants are then given by 



Since the large r behaviour of f n {^) is the same for all n (f n ~ r~ 2 as 
r — > oo), the different excited states only differ in their decay constants as 
a result of the different normalizations of the /„. 

We require that the kinetic terms for the different rho excitations are all 




4 (a 3 /;(a)) 2 



(7.2) 
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canonical which implies imposing 



I 



A 




(7.3) 



In the original AdS/QCD model (section [5.5. II) with the UV cutoff at infin- 
ity, one finds the decay constants grow as the square root of the excitation 
number n 



If one matches g§ to the perturbative high energy vector correlator [2, 113], 
giving g\ = 12tt 2 /N c then the p meson has a decay constant Fp/ 2 = 260 
MeV compared to the physical value of 345 MeV. 

In figure 17.21 we display the results of the same computation with a UV 

1/2 

cutoff present. For low cutoffs F po rises: with A = 194 MeV F%* = 478 
MeV. Comparison with the physical value again hints that a low cutoff is 
appropriate. 

On the other hand, as the cutoff is brought down the \fn behaviour 
(argued for in [147]) is apparently lost and the higher resonance decay con- 
stants fall relative to the n = case. The reason for this is that the cutoff 
impedes on the values of r where the wave functions of the eigenstates are 
substantial. By the time that the cutoff is of order a few hundred MeV the 
integral for the normalization is dominated around the cutoff. This makes 
the computation of the decay constant suspect - formally one needs a de- 
scription of the physics to higher energies which may lie beyond the region 
of perturbative validity for the supergravity. 

This contrasts with the computation of the masses in the previous section 
- those values are determined by requiring regular solutions in the infrared 
away from the cutoff. 
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A=1940 MeV 





1 12 

Figure 7.2: The decay constant F p for different p excitations, plotted for 
varying UV cutoffs 

7.4 Anomalous Dimensions at the Cutoff Scale 



There remains the possibility that the classical dimension of any operator is 
not respected at the quantum level. Here we investigate possible effects to 
the operator g7 M dj. Classically is has dimension three, and this is encoded 
in the gravity dual by fixing the UV boundary condition as 



f n (r)~r 2 , f' n (r) ~ -2r 



(7.5) 



In general though we should write 



fn(r) 



-lot 



(7.6) 



and numerically investigate the preferred value of to. The result is that for 
all values of the cutoff, from oo down to our tuned value of 194 MeV, the 
preferred value of u is two, so it appears that the classical value is conserved 
at the quantum level. 

However, this is a rather naive analysis: if one evaluates the derivative 
of the flows f n { r ) at finite values of r in the case of the model with an 
infinite cutoff, the derivatives differ for each value of n. Hence we should be 



*Thcre is a theorem [148] that says that conserved currents are forbidden from ac- 
quiring anomalous dimensions. However this theorem is only true for theories that do 
not contain vector mesons with the same quantum numbers as the conserved current, 
hence the theorem is not applicable here 
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considering the possibility of a different value of o> for each f n . Pragmati- 
cally, tuning each radial excitation is both time consuming, and means that 
we are introducing so many free parameters that we have lost all powers 
of prediction. However it is interesting to see how much the anomalous 
dimensions may contribute. The results for our best fit of A = 194 MeV is 
shown in table 17. lT a). All corrections are less than 10%, which in AdS-QCD 
terms, is negligible. 

7.5 Coupling of Higher Dimension Opera- 
tors 

The second cause for concern when introducing a finite UV cutoff in the 
strong coupling regime is the contribution of irrelevant operators. As we 
pointed out in section 17. 1[ there are an infinite number of such operators, 
and it is not known how to encode the majority into the gravity dual. 
However, one might hope the couplings of lower dimension operators would 
grow fastest as one moved into the non-perturbative regime. We know how 
to encode one simply [144,149-151], so we investigate its effect. 
Let us rewrite the metric ( 15.1 ft of model [2] as 

ds 1 = H 1/2 dr 2 + H' 1/ \ u dx"dx ,/ (7.7) 

with H equal to r -4 in 15. 11 We can deform the AdS space by allowing it to 
return to flat space asymptotically: 

H(r) -> r" 4 + a = r" 4 (l + or" 4 ) (7.8) 

The parameter a is a symmetry singlet and has energy dimension minus 
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Table 7.1: Preferred values for: a) the quantum dimension u> of the operator 
q^q (left), and b) the field dual to G Tr F 4 (right) 



four. Therefore it is identified with the coupling G of the operator Tr F 4 . 



We repeated the fit to the lightest five p mesons masses in this deformed 
geometry. The results are displayed in table 17.11 (b). 

The results show two effects. Firstly that a grows as A decreases. This 
is to be expected since a changes the large r part of the metric. To change 
the results when only a small r slice of the metric is present needs a large 
a. Secondly, the result at A = 194 MeV is an exception to this rule. In fact 
changing a in this case only has an effect on the third significant figure in 
the error. This reflects just how good the fit is from a simple cutoff. We 
conclude that at an appropriately low cutoff, the influence of the G Tr 4 is a 
small effect. 



7.6 Discussion 

A perturbative gravitational dual of QCD should only be expected to work 
at energies below a few GeV at best, where QCD is non-perturbative. We 
have investigated imposing a UV cutoff on an AdS/QCD model of the p 
mesons and found that the data has a fit at the 2% level with a UV cutoff of 
a few hundred MeV (compared to a fit of 21% with an infinite cutoff). This 
leads us to the conclusion that the holographic description of QCD should 
only be used at low energies on a quite small radial interval. However this 
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conclusion is countered by table l7TIT a): all the errors are small, and so it is 
also plausible that the holographic description of QCD could stretch upto 
several GeV. This would dovetail nicely with chiral perturbation theory and 
perturbative QCD. 

We have also looked at fitting corrections to the anomalous dimension 
of the operator q^^q and introducing a coupling of the operator Tr F A . 
Although these corrections could be used to fine tune the fit by a percent 
or so they do not appear to be significant corrections to the model. Of 
course these are only easily implementable examples from an infinite set of 
possible corrections but finding the corrections to be small provides further 
understanding of the success of the basic AdS/QCD models. One could 
also try to include the vacuum expectation values of more operators in 
the metric (see for example [152]) and a dynamical, predictive mechanism 
of chiral symmetry breaking [153]. Such effects would be important to 
study the pion and axial vector meson sectors of the model. As explained 
in [117] the model used here does not give a good prediction of these sectors 
because the dilaton form, put in to give the \fn rise in masses, does not 
lead to a sensible condensate prediction. If one attempted to tackle all of 
these problems then most likely the number of free parameters would rise 
faster than the number of available data points. Of course this reflects the 
fact that a perfect action is in the end just a reparameterization of the full 
QCD spectrum. We hope though that we have identified the imposition 
of a UV barrier as an important correction and that these other effects are 
sub-leading in the p sector. Putting together a complete model of all sectors 
including the baryons remains as an important challenge. 



105 



Chapter 8 



Hadronization 



We now turn to a notoriously difficult problem of strongly coupled QCD: 
hadronization. Hadronization is the process of the formation of hadrons 
from gluons and quarks. It is widely believed that this is how hadrons were 
formed in the very early stages of the universe, as the quark-gluon plasma 
cooled sufficiently. Nowadays, those kind of energies and densities are only 
reached in extreme natural conditions, such as in the cores of neutron stars. 
They can also be reproduced in the laboratory, typically at high energy 
particle colliders of the type found at CERN, FermiLab or DESY, and most 
recently with the Relativistic Heavy Ion Collider (RHIC) in Brookhaven, 
USA. 

All these colliders are based on the simple premise of colliding two par- 
ticles together at very high energies, and seeing what the result is. Richard 
Feynmann likened the process to smashing two Swiss watches together to 
figure out how the watches work. It is a crude method, but has been very 
successful. Of primary interest are the hard processes which occur immedi- 
ately after the collision: these processes can be calculated using Feynman 
diagrams and have been used to verify the predictions of QCD to within 
an error of a few percent. However, the vast majority of the collisions after 
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the initial hard collision are soft, and cannot be treated using perturbative 
QCD. The question we address in the this section is one of the most simple: 
"If we annihilate an e + e~ pair at high energy, what particles, and in what 
numbers do we detect in our apparatus?" This simple question hides a 
multitude of very difficult issues. If taken at its most fundamental interpre- 
tation, it would require us to predict the mass spectrum of all hadrons, their 
branching ratios, as well as the physics through a probable phase change. 

Instead we limit our answer to just predicting the initial yield of hadrons 
following the (assumed) annihilation of the e + e~ pair. For the hadronic mass 
spectrum we use what has been phenomenologically observed [24], and we 
use the same source to obtain the branching ratios of every hadron. Our 
treatment is not without precedent. [3] assumed that after the quarks freeze 
into hadrons, they may be described as a hadron gas in thermodynami- 
cal equilibrium. [3] made no attempt to predict branching ratios or mass 
spectra. 

It will be prudent to first review the elementary facts of hadronization. 

8.1 Basics 

The essence of hadronization is shown in figure 18. 1[ with the black arrow 
in the middle symbolising the little understood hadronization process. A 
computer re-enactment of a single e + e~ collision is shown in figure 18.21 

It shows that for any one collision, many hadrons are produced. The 
precise number of hadrons produced in each collision varies, reflecting the 
probabilistic nature of quantum field theory. But over many collisions, we 
can compile an average number per collision for each hadron. A sample of 
these results is shown in table IHTT1 (we tabulate the complete results later). 
The heaviest hadron detected, the Q baryon, has a multiplicity of 0.0014, 
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hadrons (ti ± , K ± ,p,A° etc.) 



Figure 8.1: Hadronization 



meaning that for every 10,000 collisions, it is seen fourteen times. Over the 
same number of collisions, we would have seen 918,000 7r°s. 

The creation of the qq pair in figure 18.11 immediately after the e + e~ 
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Figure 8.2: A re-enactment of a single e + e collision detected at CERN 
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Table 8.1: Mean multiplicities for a typical e + e collision at y/s = 91.2 
GeV [3] 



collision is a hard process, and is easily modelled by Feynman diagrams 
such as that shown. What as yet is unknown is how to model the process 
represented by the black horizontal arrow of figure 18.11 Although there 
may still be some hard processes occurring, the vast majority will be at low 
energies, especially as the quarks form more and more hadrons. These soft 
processes represent a difficult conceptual challenge. 

Conservation of momentum means that the qq pair will be travelling 
back to back. Conservation of energy means that they will be travelling 
very fast. Colour confinement forbids the quarks from existing individually, 
and how they form colour singlet composites is the focus of several theories. 



In the independent model the qq act independently and each com- 
bines with quarks and antiquarks spontaneously created from the vacuum 
to form hadrons. 



Alternatively, the Lund string theory treats all but the highest-energy 
gluons as field lines, which are attracted to each other due to the gluon self- 
interaction and so form a narrow tube (or string) of strong colour field. This 
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colour tube forms between the qq pair, and fragmentation of the strings in 
the tubes causes the formation of hadrons (see [154] for more details). The 
Lund theory is part of a more general class of parton fragmentation models 
used by event generators such as PYTHIA and JETSET. However these are 
complex multi-parameter Monte Carlo event generators that do not reflect 
the simplicity of QCD, and have little predictive power. 

A recent model which has had striking success is that of [3]. It assumes 
that each quark/ ant iquark from the qq pair forms a hadronic gas in thermo- 
dynamical equilibrium. The gas then cools and expands, until it is no longer 
in equilibrium, at which point the hadrons are frozen out, and these are then 
the particles that travel to the detector. So this model only requires three 
parameters: the temperature and volume of the gas at the point of freeze- 
out, and a strangeness suppression factor to account for the superior mass 
of the strange quark. The strangeness supression factor should, in principle, 
be derivable from the known mass of the strange quark. However, such an 
attempt is not made in [3], and 7 S is left as a free-fit parameter. (Factors for 
the charmed and bottom quarks can also be introduced, but these quarks 
are so massive, their production rate is negligible). This model can repro- 
duce the results of LEP and PEP-PETRA to 34% and 38% respectively. 

All models of hadronization have two parts: predicting the initial yield 
of hadrons directly after annihilation, and then allowing for decays of those 
particles in transit towards the detector. Both of these processes are QCD 
processes, and so should be fully predictable from the theory. However 
just predicting one of these two processes is a major challenge: here we 
concentrate our efforts on the first stage. We will then use experimentally 
measured values for the branching ratios of all the hadrons [24] to predict 
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what would actually be seen in a detector. 

8.2 The premise 

All AdS-QCD models provide a weakly coupled five dimensional gravita- 
tional theory that describes the hadrons of QCD. Each hadron and its ex- 
cited states (e.g. the 'stack' of p meson masses: p(770), p(1450), p(1700) . . .) 
is described by a five dimensional field that shares the Lorentz and global 
symmetries of the hadron. In the gravitational theory one seeks solutions for 
those fields that separate the 3+1 dimensional dependence from the extra 
radial dependence. More succinctly, we seek solutions of the form g n {r)e lknX 
with = M%. There are only regular solutions on the space for discrete 
masses: these correspond to the meson masses. The functions g n {r) form 
an orthogonal basis with the appropriate weighting function. 

We will assume that at the point where the quarks form hadrons, the 
energy of the event is democratically available to all hadronic channels. 
We describe the initial condition as some deposition of energy into the 5D 
model's stress-energy tensor. The radial dependence of the stress-energy 
dependence can be expanded in terms of the functions g n and will deter- 
mine the relative multiplicities of each particle in a hadron stack. The x 
dependence will determine the energy and momentum of the hadrons. In 
this paper we will simply concentrate on the multiplicities. 

The result for the multiplicities depends on the choice of the function 
expanded in terms of the g n and this represents the matching to the under- 
lying asymptotically free QCD dynamics. The simple guess we will employ 
is that we should treat all hadronic channels equally and pick a Gaussian for 
the initial condition. The height of the Gaussian determines the absolute 
value of each particle's multiplicity and hence is a free parameter which we 
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fit (this parameter can be re-expressed as the average energy per hadron, 
k, which we detail later). The width is also a free parameter which we fit. 
Finally the data displays a suppression factor on the production of strange 
quarks - this is not surprising since they are more massive. We reflect this 
factor in the underlying dynamics by a strangeness suppression factor 7 S 
which multiplies the Gaussian for each strange quark in the hadron (note 
that if mixing occurs the strange quark content need not be an integer). j s , 
which is also used in the thermal hadronization models [3], is another fit 
parameter in our analysis. 

For any collision with an initial energy of y/s, no particle with a rest mass 
greater than \/s will be produced, so in principle we would need a cutoff 
which would depend on the centre of mass energy of the initial collision. In 
practice we shall simply include all known states with mass below 1.7 GeV. 
Above this value the experimental data on the full spectrum becomes patchy. 
In addition, the high mass states in this range are only produced with very 
small multiplicities and have a minimal effect on the lighter particle results. 
Since the do (980) triplet is widely thought to be a bound kaon state [24] 
we include this as a decay channel, but not as a state that can be formed 
immediately after the annihilation. 

The final ingredient we require is a specific AdS/QCD model of the QCD 
hadrons. We will adapt a string theory derived model of chiral symmetry 
breaking, based on the Constable- Myers geometry of chapter H] that includes 
the vector mesons and pions [69]. That both the vector mesons and the 
pions are included is an important feature: the mass spectra of the pseudo- 
Goldstone bosons are significantly different to all other hadrons in QCD. We 
will then assume that the g n functions associated with each hadron stack are 
not that different from the p-stack functions - we will simply reproduce them 
but with the mass of the lightest stack member tuned to the experimental 
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value. Similarly the pion stack can be used to reproduce the towers of states 
associated with each pseudo-Goldstone of the chiral symmetry breaking (i.e. 
the pions, kaons and eta meson). The relative weighting is parameterised 
by R, a measure of gyuN, which we fit. 

8.3 Holographic Hadron Basis Functions 

Our holographic model of the hadron spectrum is based on the Constable- 
Myers geometry discussed in chapter HI We then compactify it to five di- 
mensions, in exactly the same way we did in chapter [6] when trying to 
regulate the infrared of holographic descriptions of QCD. We remind the 
reader that this model spontaneously forms a chiral condensate, and has a 
complete orthogonal set of functions for both the vector meson sector and 
pionic-like sector. Henceforth we denote these two sets as g n (r) and f n { r ) 
respectively, matching the notation used in chapters 0] and El 

The mass of the lowest lying p meson can be dialled by choosing the 
conformal symmetry breaking scale in the model, and the mass of the lowest 
lying pion can be dialled by choosing the asymptotic quark mass. The 
excited states in both stacks are then predicted: m p * = 1737 MeV, m^* = 
1701 MeV (c.f. experimental values of 1459 and 1300 MeV respectively). So 
whilst they don't precisely reproduce the experimental values the pattern 
is at least roughly right. 

As stated above, we assume that the g n functions associated with each 
hadron stack are not that different from the p-stack functions, and simply 
rescale the r coordinate such that the mass of the lowest member of each 
stack is correct. 

We class the pions, kaons and 7/(548) mesons as pionic, and use the f n 
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Figure 8.3: The normalised holographic hadron basis functions for the 7r° 
stack (left) and the p meson stack (right) 

functions to model them. For each pionic meson, we dial the asymptotic 
quark mass such that the mass of the lowest member of each stack is correct. 
To give an idea of the functions that result, we plot the first three functions 
for f n ,g n in figure El 

It must be clear to the reader that this situation is not ideal: our ad hoc 
model predicts over massive excited states, and we have crudely broadened 
the scope of a model which originally contained only the vector meson sector 
and a single pionic sector. We would prefer a model which contained all the 
fields of QCD, and predicted their masses correctly. However, this would be 
a major breakthrough in progress towards a holographic dual of QCD. This 
breakthrough has yet to be achieved, and so in the meantime we proceed 
with the model as presented here. 



8.4 Overlap Computation of Multiplicities 

With our holographic functions f n , g n in place for each hadron stack we can 
now proceed with computing the expected initial yield in a hadronization 
event . 

We assume that all the five dimensional fields \l/(r) have a common 
initial condition of a Gaussian centred at r = and of width A. To find the 
multiplicities of each stack member we compute 
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vj/(r) w(r)g n (r)dr (8.1) 



o 



where w(r) is the weighting function associated with the basis functions g n . 

The multiplicity is simply given by (? n multiplied by (2J+1) where J is 
the spin of the hadron. 

There £1X6 cL number of special cases, which we address in turn. 

Pseudo-Goldstone Bosons 

The pseudo-Goldstones are described by a separate holographic field, 9, in 
our model and are represented by functions f n { r ) ■ The relative contributions 
the fields make to the stress-energy tensor are 

T rr ~ A 1 (r)(9 r ^) 2 + A 2 (r) V ^(d r A,)(d r A u ) + ... (8.2) 

Ai, A2 are functions of r, completely calculable from the holographic model, 
using the usual definition of the stress-energy tensor: 



T/w = ~ 2 (V - fi^mattcr) (8-3) 



To ensure the fields and 6 see the same contribution to the stress- 
energy tensor we must rescale the Gaussian. If our standard Gaussian is 
\l/(r) then for the energy to be equally divided between all sectors we require 

d r 9 = y/ftdr*. 

Strangeness Suppression Factor 

Since the underlying asymptotically free dynamics will distinguish the strange 
quark from the up and down quarks, we also multiply the Gaussian by a 
factor of (7s) CT where a is the strangeness content of the stack. 7 S is then the 
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second fit parameter in our procedure. Note that this procedure is rather 
crude because different members of a stack may mix to varying degrees with 
other states. For example the 7/(548) has 32% strangeness content while the 
r]** has 100%. In these cases we set 7 S by the strangeness content of the 
lightest member of the stack. 

Height of Gaussian 

The normalisation of the Gaussian tells us the relative multiplicities of the 
various hadrons in an event. An overall multiplicative factor k sets the 
absolute number of each species and we fit this value, k determines the 
total number of final state particles (before allowing for decays in transit to 
the detector), and hence we express it as the average hadron energy in the 
collision. 

A Fourth Parameter 

Our choice of holographic dual also contains a free parameter, R, which sets 
the 't Hooft coupling in the gravity dual. We fit it to the data. However, R 
is not in the same class as A, 7 S , k. R has a sound theoretical background, 
and would not be present if the holographic dual to QCD was known. 

Decay in transit 

Once we have calculated the initial yield of hadrons, we then have to al- 
low for decays of the particles in transit from the interaction point to the 
detector. Branching ratios are taken from [24], and particles that can be 
detected at LEP (whose results we will compare to) are set as stable. All 
the other particles are allowed to decay through the decay channels until 
they reach one of the stable particles. In this way we get a list of numbers 
which is what our model predicts would be seen at LEP. 
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8.5 Predictions 



We compare our results both to e + e collisions performed at LEP (-^/s = 
91.2 GeV), and at PEP-PETRA (y/s ~ 30 GeV). Average multiplicities of 
various hadrons have been compiled in [3], and we reproduce them here, 
along with our results, in table 18. 21 For both sets of results we have per- 
formed a four parameter fit so as to minimise the rms error. The pion 
predictions are expected to be low because our holographic model predicts 
over-massive excited states. Thus yields of particles such as 7r(1300) and 
p(1450) are unnaturally suppressed which would otherwise be expected to 
give significant contributions to the pion multiplicities. Despite this, the 
fits are very good, with rms errors of 47% and 53% for ^/s = 91.2 GeV and 
\/s ~ 30 GeV respectively. The Q and rj' yields aren't ideal, but previous 
models have had the same problem [3], especially when matching to the 
PEP-PETRA data. Our model is expected to have poor predictive power 
for these two hadrons: the 7/ is technically a pseudo-Goldstone boson, but 
instanton effects cancel out this effect. Our model just treats it as non- 
Goldstone boson, which is probably over simplistic. In addition, it mixes 
heavily with r/(548): our model contains no good parameterisation of mix- 
ing. The Q baryon also highlights weaknesses of our model: it is a baryon, 
not a meson (our functions g n originally source from a model of mesons, 
not baryons); it is very heavy, and so is strongly affected by our exclusion 
of hadrons above 1.7 GeV; it contains three strange quarks, making it very 
sensitive to j s . 
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Table 8.2: Results of the model for hadron yields at yfs = 91.2 GeV (centre 
column) and y/s ~ 30 GeV (right column). The relevant 3+1 free parameter 
values are shown in the final row. 
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0.0089 


0.0083 




0.012 


0.012 


0.0088 


0.0083 




0.0040 


0.0033 


n/a 


n/a 


n 


0.0011 


0.0014 


0.0008 


0.007 


A(MeV),K(GeV),7 s , J R 


150, 4.96, 0.97, 2.6 


152, 2.35, 


0.97, 2.4 



8.6 Conclusions 

We assumed that every hadron in QCD can in principal be represented by 
a function in the r coordinate of the 5D holographic theory of QCD. We 
then proposed that hadronization can be modelled by hypothesising that 
the initial yield (that is before the particle created starts decaying) for any 
hadron is given by the square of the overlap between the function which 
represents the hadron, and a Gaussian, centred at the origin, with a width 
of A. In addition we have two other parameters in the theory; a strangeness 
suppression factor to account for the heaviness of the strange quark, and k 
which determines with what energy the particles leave the interaction point. 

With the full holographic dual to QCD currently unknown, we made 
some ad hoc assumptions to achieve a full set of functions which represented 
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every hadron. We could then put our hypothesis to the test, and compared 
the results to e + e~ collisions made at LEP and PEP-PETRA. The results 
were surprisingly good. 

It is also interesting to note that although we included a strangeness 
supression factor, 7 S , in the manner of [3], we found it to be almost equal 
to one in both fits. This suggests that in practice such a parameter is not 
necessary in our model. In addition, our holographic model of QCD predicts 
over-massive excited states, and this probably leads to a suppression of the 
pion yields in table [821 A holographic model with a more accurate hadronic 
mass spectrum would almost certainly give better results. 

The model of hadronization presented here is applicable to all particle- 
antiparticle annihilation events, where the fireball after the collision has no 
residual quantum numbers. Broadening this model to include events such 
as deep inelastic proton-proton scattering, and heavy ion collisions would 
clearly be desirable. We leave such an analysis for the future. 
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Chapter 9 



Conclusions 



QCD is far from being fully understood. The non-perturbative aspects of 
the theory are particularly intriguing: confinement, asymptotic freedom, 
chiral symmetry breaking, the hadron spectrum, hadronization, and the 
QCD phase diagram all offer exciting research opportunities. The AdS / CFT 
correspondence has proved to be a valuable tool in our understanding of 
strongly-coupled gauge theories, and the related AdS/QCD models have 
reproduced many aspects of QCD very successfully. 

In chapter 1 we described some of the more beguiling aspects of QCD, 
and reviewed the very basic fundamentals of string theory. In chapter 2 
we introduced the AdS/CFT correspondence, in particular explaining how 
every operator in the field theory can be uniquely matched with a field on 
the gravity side. We used chapter 3 to show how to include flavour into the 
correspondence, via a D7 brane probe. This gave a relationship between the 
quark mass and quark condensate, and we showed that in this geometry, 
chiral symmetry breaking could only be induced by massive quarks: the 
vacuum does not spontaneously break the chiral symmetry. We went on to 
calculate the (scalar) meson spectrum in this set-up by looking at the scalar 
fields on the surface of the D7 brane. We also showed that the vector meson 
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spectrum in this geometry was identical to the scalar meson spectrum, which 
was to be expected by supersymmetry. 

In chapter 4 we introduced a geometry, called the Constable-Myers ge- 
ometry, which was singular and non-supersymmetric. It is an important ge- 
ometry because it was the first known geometry to holographically demon- 
strate chiral symmetry breaking. We summarised its salient properties: 
confinement and the demonstration of a mass gap (via the glueball spec- 
trum). We then introduced quarks with a D7 probe, and confirmed the 
spontaneous breaking of chiral symmetry. The pion and vector meson spec- 
tra were calculated, and we plotted m p against m^. This was compared 
against equivalent lattice data for QCD at large N, and the results were 
found to be in excellent agreement. 

Chapter 5 brought us into the realm of bottom-up theories, generically 
called AdS-QCD models. We studied in detail the most famous example of 
such a model because it was historically important, simple to follow, and 
very successful. By including just three operators in a truncated AdS space, 
it was possible to holographically model chiral symmetry breaking, and 
hence predict the masses and decay constants of the a\ meson, the p meson, 
and the pion. The predictions were to within about 10% of the true QCD 
values. We then generalised our first AdS-QCD model to include a dilaton in 
the action. This had the welcome effect of making the vector mesons follow 
the Regge trajectory of ~ n (c.f. m n ~ n without the dilaton), but the 
unwelcome effect of ruining the spontaneous chiral symmetry breaking of 
the model, and thus removing any meaningful predictions of the axial vector 
spectrum. We then briefly described the many other phenomena which have 
been described in AdS-QCD models. 

In chapter 6 we attempted to address some of the shortcomings of the 
existing AdS-QCD models: the lack of asymptotic freedom, the hard IR 
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cutoff, and no dynamical determination of the quark mass and condensate. 
We took the Constable-Myers geometry of chapter 4 and truncated it by 
removing the S 5 part of the spacetime. This then provided a model which 
was everywhere smooth, and had all the features that we analysed in chapter 
4: confinement, a mass gap, and spontaneous breaking of chiral symmetry. 
The model, like QCD, only depended on two parameters: the light quark 
mass m q , and the strong coupling scale A (c.f. the model of chapter 5 which 
had three parameters). For comparison with chapter 5, we calculated the 
same quantities and found agreement with experiment to the 10% level. 

Whilst chapter 6 had regularized the IR of AdS-QCD models, chapter 
7 dealt with regularizing the UV of such models. Without regularization in 
the UV, we were implicitly trying to solve a string theory in the strongly 
coupled regime. The simplest approach was to include a hard UV cutoff 
in the gravity dual. This removed the strongly coupled regime from the 
theory, and made predictions of the vector meson spectrum to within 2% 
of accepted values. However, the UV cutoff complicated the operator-field 
matching process. Firstly if a field had an anomalous dimension it would 
have to be fitted to a different operator with the relevant dimension, and 
secondly irrelevant operators would become more and more dominant as we 
neared the non-perturbative regime. Hence they would have to be encoded 
in the gravity dual. Without being able to address each complication fully 
we instead looked at an example of each and found that, for these particular 
examples, the effect was not large. 

The final chapter addressed hadronization. We postulated that after 
some annihilation event in a particle accelerator, the residual energy took 
on the shape of a Gaussian. The initial multiplicities of some hadron species 
would then be given by the overlap between this Gaussian and the radial 
function which represented that particular hadron in the QCD holographic 
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dual. Once the initial multiplicities were known they were then allowed to 
decay according to their published branching ratios until they reached the 
detector. With the full holographic dual to QCD currently unknown, we 
used the model of chapter 6. We compared our predictions to e + e~ collisions 
at centre-of-mass energies of 30 GeV and 91.2 GeV. Data which ranged over 
four orders of magnitude were predicted to the 50% level. 

The overall conclusion to be drawn from this thesis is that the AdS- 
QCD correspondence is in its infancy, but can still reproduce many aspects 
of QCD at very surprisingly accurate levels. Regulating the ultraviolet and 
infrared regimes of the dual are necessary steps that don't disrupt the good 
predictions. Most excitingly, it seems that AdS-QCD models can model 
hadronization events very well with only a few free parameters. 
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